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1. Introduction 

In the preceding paper ^ (referred as I), we have developed a general formalism 
based on Green functions to calculate the electromagnetic field scattered by a random 
medium with rough boundaries. In using these Green functions, we have determined 
the average electric field, also named the coherent field. In this paper, we use the 
formalism developed in I to calculate the diffusive intensity also named incoherent 
intensity. 

The study of electromagnetic wave propagation through random media has been 
an intensive field of research for many decades [21131111311011313011311111111131 

1141 1151 . Many aspects of the transport of waves in such media are well described 
by the phenomenological radiative transfer theory [3131311113 ^1 ^3 ^] ■ In this 
approach, the main quantity, used to describe the propagation, is the specific intensity 
X(i?, k) (if we take into account the polarization of the wave, the specific intensity 
can be defined as a Stokes vector or a tensor) which gives the power flux per unit area 
and solid angle at the point R which goes in the direction k. In writing a balance 
equation on the energy, it is shown that the specific intensity satisfies a Boltzmann 
type equation called radiative transfer equation. If the particles are inside a slab with 
a permittivity different from the surrounding medium, boundary conditions must be 
added to the radiative transfer equation in order to calculate the specific intensity. For 
rough surfaces, these boundary conditions are expressed with scattering operators, 
where several approximate analytical expressions exist depending on the roughness 
of the surface. Numerical calculations of the radiative transfer equation taking into 
account rough boundaries can be found in references [3 El El ITHj 

However, all previous studies on the radiative transfer theory are based on 
heuristic principles since in these works, the specific intensity is a fundamental quantity 
which is not defined from the electromagnetic field E and B. The link between the 
electromagnetism and the classical radiometry theory is due to Walther ^Hj who has 
first recognized that the specific intensity can be deduced from the Wigner function 
I{R,k) ^ Const X J d^r exp^-ik ■ r) E{R + r/2) ^g) E* (R- r /2), where the tensorial 
product (g) permit to take into account the different polarizations of the waves. If 
we suppose that the electric field depends slightly on the coordinates R compared 
to r (quasi- uniform fields hypothesis ^S])) it is shown that the Wigner function is 
written as I{R,k) = Const x (5(||fc|| — Ko)X{R,k), where k = fc/||fc|| and Kq the 
wavenumber. In this case, the function I{R, k) can be identified with the specific 
intensity. Nevertheless, there are some differences between this definition of the specific 
intensity with the radiometric one. In fact, this electromagnetic definition of the 
specific intensity is not always a positive function, which is in contradiction with the 
radiometric interpretation in terms of power. This problem has been the subject of 
much study J,3, 20, 21 where it has been demonstrated that in the geometrical limit, 
the electromagnetic definition of the specific intensity is always positive. In other 
cases, the negative value of the specific intensity are due to interference effects. 

As soon as we have determined the relationship between the electromagnetic 
field and the specific intensity, we can derive the radiative transfer equation from 
the Maxwell equations. The procedure is to write the Maxwell equations in an 
integral form with the help of Green functions and to apply the Wigner transform 
to the equation derived. Then, in differentiating this equation we obtain the radiative 
transfer equation [2a[23[21[23[23[2a[23l23l2IlinilSai33ISIISS|- Furthermore, in 
starting from the wave equations, we are able to take into account new contributions 
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to the scattered intensity such as the enhanced backscattering and the correlations 
between the scatterers that can not be introduced in the phenomenological radiometric 
approach. The objective of this paper is to derive the radiative transfer equation 
from the wave equation by taking into account rough boundaries. Several works 
have investigated this topics, but they describe the scattering by the rough surfaces 
either by using the small-perturbation method |3()l 1371 138| or in an unconventional 
fashion [211 Ell • our approach, we use scattering operators which are a versatile 
and unified way to describe how the wave interacts with the boundaries 01]. To 
use these operators, we have introduced two kinds of Green functions (I). The first 
one Gsv{''','''o) describes the field scattered by the volume (V), which contains the 
scatterers, and by the rough surfaces (S). The second type of Green function is 
Gs{r, ro), which describes the field scattered by a slab with rough boundaries where 
the scatterers have been replaced by an homogeneous medium of permittivity ee named 
effective permittivity. As demonstrated in I, the Green functions Gs{r, Tq) are easily 
expressed as a function of the rough surface scattering operators. In using these 
kind of Green functions, we were able to separate the contribution coming from the 
surface and the volume. The main advantage of our approach is that the equations 
obtained are similar to the equations generally used to describe the wave scattered by 
an infinite random medium 9, 24, 25, 27| l^ . To take into account the boundaries, we 
replace the Green function of an infinite homogeneous random medium with Gs{r, Vq). 
Yet, there is a slightly difference from the classical procedure. Usually, we use the 
normalized vector k (||fc|| = 1) to describe the propagation direction of the wave since 
the generalized Dirac function in the specific intensity X(r, k) imposes that k = Kq k. 
In this study, we will use the two-dimensional vector p = k^e^ + kykz to describe the 
propagation direction, and we will recover the vector k by following decomposition: 
K^k = p + aao{p)ez, where ao{p) ~ {K^ — p^)^/^, and a is the sign of the vertical 
component of A; (a = sgn(A; • 62)). This choice, which is unusual for the radiative 
transfer theory, is in fact the standard in scattering by rough boundaries where the 
vertical axis z has to be differentiated from the axis x — xe^ + ye.y for a surface 
profile defined by z = h{x) |41) . Furthermore, the distinction between upward waves 
for a = + and downward waves for a = — is useful to write the boundary conditions. 
At the end of this paper, we will explain how to rewrite the equations obtained in the 
usual form with k vectors. 

2. Cross-section 

The geometry of the problem and the notation are described in paper I. In order to 
characterize the scattered intensity by an object, we usually introduce the bistatic 
cross-section, which is the power scattered per solid angle normalized by the incident 
power flux. In this paper, we will use a generalization of this concept called MuUer 
bistatic cross-section, which permits an accounting for every state of polarization of 
the incident and scattered waves. First, as an intermediate of calculation, we have to 
introduce the scattering operators describing the field scattered by the rough surfaces 
and the random media ||. For an incident plane wave, 

E°'{r) = E°\pq) e^Po-^'io^oiPo)^ ^ (1) 

II Notice that in the paper I, we have introduce the scattering operators describing the field scattered 
by an homogeneous slab with rough boundaries. However, we can also define, in the same way, 
scattering operators when the medium is inhomogeneous. 
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the field scattered by the random medium and the rough boundaries is: 



where ao{p) = y/ Kq — p^, and using the notation defined in I, we decompose the 
vector £J°*(pq) and the dyad Rsvip\Po) on the following basis: 

E°\p,) J2 E'>\po)0el-ip,) (3) 

Rsvip\Po)^ J2 ^sv(p|Po)/3,/3oe|^"^(p)47(Po): (4) 

where the polarization vector e|^^(p) for the polarization TM (/? = V) and the 
polarization TE (/? = H) are defined by 

e"/(p)=±^p-^e., e,ip)=e,.p. (5) 

Far from the scattering medium (i^Tol 1^1 1 3> 1), we can obtain an asymptotic expression 
of the integral in equation |5J) by using the stationary phase approximation p2j : 

piKo\\r\\ _ 

Efvir) = -j^ f{p\Po) ■ E"\Po) , (6) 

with 

P = ^o^, (7) 

7{p\Po)-^Rsvip\Po)- (8) 
From expression we can decompose the vector E'^\^{r) on the following basis: 

E%ir)^E%ir)ve'+ip)+E%ir)HeHip), (9) 

where the vector p is defined by {Tj). The incident and scattered intensity can be 
described with the help of Jones tensors (13) : 

J°%r) ^eoE%{r)^E'^^^-{r), (10) 
= eo Yl Efy{r)pEl\;{r)prel+{p)®elt{p), (11) 

P,13'=H,V 

J^\p,) = e,E\p,)®E'"{p,) (12) 
= eo i^"(Po)/3i?°"(Poke°;(po)®e°7(po), (13) 

often written in matrix form I43| : 

^^)\-^^yEfy{r)y[E%{r)HY \E%{r)y\^ j' ^^^^ 

[J y^'n-^^ y E^\r)y[E''\r)H]* \E°'{r)y\'' J' 

Here 

eo = 7, (16) 
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and e^acy c.^ac are, respectively, the permittivity and the speed of hght in the vacuum, 
and no = ^/eo the optical index of the medium 0. The incident and scattered Jones 
tensors are related by the tensorial cross-section defined by 

4-7r 

nP\Po)--J"\P,)^^\\r\?J''{r). (17) 

where A is the area illuminated by the incident wave and : is the product between two 
tensors as defined in [Appendix A| In the vectorial basis [e^{p )e/f (p)], we have 

A'jr 

E ^{p\p,)0p-A,P',El{p,)Ef^{p,) = ^\\r\\'E%{r)pE%*{r)p, (18) 

l3o,l3'o=H,V 

where 

^{p\Po) = <^{p\Po)pp';PoP'o{f'l^{p) ® e°^(p))(el^;(Po) ® e°7(po)) • (19) 

!3,l3'-l3o,l3'o=H,V 

Definition is obviously a generalization of the usual scattering cross-sections since 
the elements a'vv-yv, '^vV;HH, o'hh-.vv, and 'Whh-.hh of the tensor '^pp'-^p^p'^ are 
given by 

47r \\r\?\Efy{r)v\' 
"^^^^^ - T \E^-^{p,)vV ' ^^^^ 



47r \\rmE%{r) 
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47r \\rW'\Efy[r)H? 

"^^^^^ - T \e^Kp.)v? ' ^ ^ 

4^ \\r\?\Efy{r)H? 
^ T \E-^{Po)H? ' ^''^ 
which are the usual scattering cross-sections per unit area cyy, cryn, crij\/ and ffuH- 
From equation 0, we deduce that 

^"'(^) = i(7(p|Po) ® 7*(P|P0)) : .7°XPo) (24) 



4-7r * 

5^(p|Po) = ^/(p|Po)®/ (PlPo), (25) 



if the medium is non- absorbing (eq G K+). The cross-section cr is then 

Itt 
A 

= ^Rsy{p\p,)®%y{p\p,) . (26) 

For a random medium and rough surfaces described statistically, we usually separate 
the scattering contribution into a coherent W'^°^ and an incoherent part a"''^"^: 

-^'°\p\Po) = 4^ « Rsv{p\Po) ^sv ® « \Rsv{p\Po)] »sy , (27) 

A TT 

5^™-''(p|Po) = 4^ \ « fisy(plPo) ® Ssv'(pIPo) »sv 



< fisv(p|Po) >sv < 'Rsv{p\Pg) ->sv 



(28) 



where ^^sy is the average over the surface and the volume disorder. For statistically 
homogeneous random medium and surfaces, the average of the scattering operators 



—tncoh 

we introduce the tensor Tc^^qv^ such that 
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Rsv{p\Pq) contains a Dhac distribution ( [Appendix B| ), and we define a tensor 
-R«;sv> by 

(27r)2^(p-Po)H«sy»(Po) =<^sv(p|Po) >sy • (29) 
Similarly, the average of the tensorial product also contains a Dirac distribution, and 

icoh 

{2nyS{0)n'^^^sv»{p\Po) - < Rsv{p\Po) ^Rsv{p\Po) >sy 

- < Rsvip\Po) >sv «) < fisy(p|Po) >sy • (30) 
Accordingly, the coherent and incoherent bistatic cross-section are 

W^^^iplPo) = ^ {2n)^6{p - Po)R«sv»{Po)<^R«sv»{Po) , (31) 

TT 

^ (PlPo) '^(»«:sy»(PlPo) ' (32) 

TT 

where we have used the fact that 5(0) = y4/(27r)^ for a finite patch of area A 
( [Appendix B| |. The coherent component 5"^°'' is directed only in the specular direction 
due to the Dirac distribution. In paper I, we have described how to calculate the 
coherent electric field, and we have obtain that 

« E^'iyir) S'^'XPo) • -E"^(Po) e'Po-+-°(Po)^ (33) 

Comparing equation (|33|) with the average of equation ((21 and by using definition 
(|29|l . we obtain 

S«sv»(Po)-^™\po)- (34) 
and the coherent cross-section is 

a^°\p\p,) = *1 (27r)2<5(p - p,W°\p^) ® S'°\p^) . (35) 

TT 

As was demonstrated in I, the dyad S'^°'\pq) is related to the average scattering 
matrix S {p\Pq) by the following relationship: 

< S°+°1p|po) >5- (27r)2<5(p-po)r°'(p„). (36) 

The dyad s'^^'^ {p\Pq) describes the fields scattered by an homogeneous slab of 
permittivity eg with rough boundaries. The Dirac distribution in equation Ij3fi(l is 
due to the statistical homogeneity of the rough surfaces. The coherent component 

a'^"'^ is thus totally determined by the dyad (pg) and the effective medium Ce 
described in I. Accordingly, in the rest of this paper, we will focus only on the 
calculation of the incoherent part of the scattering. As mentioned in the introduction, 
the fundamental quantity in the radiative transfer equation is the specific intensity. 
We define X™^°'*(i2, fc) as the Wigner transform of the incoherent intensity, in the 
medium 0, scattered by the random medium and the surfaces: 



k) = dVe-''=- [ « E'sUR + |) ® E'sViR |) »sv 



«£;^t'(^+|)»sy®«-E^y*(«-^)»si^ ■ (37) 
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In introducing definition (j^J in equation (|37() . and by using equation ()32|l with the 
quasi- uniform field approximation we obtain: 

^^r"-''(i?,fe) = (27r)J(fc, -ao(p))5^^"^°'^(p|Po) :^°'(Po), (38) 

TT 

with k = p + e-z- We have supposed that the medium is non-absorbing, and 
thus, the intensity I™'^°^{R,k) does not depend on R as it appears in equation 
(|2Hl- From the result in we can determine the bistatic cross-section 

by calculating the specific incoherent intensity defined by the Wigner transform H37() . 
The Dirac distribution in equation 1)38(1 insures that the vector k has a fixed normed 
||A;|| = Kq. The tensor I™'^°^{R,k) is not exactly the specific intensity used in the 
radiometry theory since it is not homogeneous to a power per unit of area and solid 
angle. The usual specific intensity T™^°''(i2, fc) can be deduce from (i?, k) from 

the following relationship: 



r"-''(i?,fc) - '—^5{\\k\\ - K^)T'-'^°\R,k) , (39) 

^0 



where k — fe/||fc||. The Dirac distribution i5(||fc|| — Kq) insure that the wave 
propagation direction fc has a fixed norm given by Kq. As demonstrated by equation 
H38|l . this property can also be written under the following form: 

T^-°^{R,k) = (27r)(5(fc, - ao{p))r^''"'{R,p) , (40) 

where 

T"^"\R,p) = -^w^-^^^pIp,) : jO''(po) . (41) 

It can be easily checked that we have the following relationship between X™™''(ii, fc), 
X'"'=°''(i?,fc) and I""'"'' {R,p): 

'^^'^■r"™''(i2,fc) = / d^fc J"™''(H,fc) = / -^X'"=°'«(H,p), (42) 



(27r)3 ' ' ' J ' ' J (27r) 

where d^fc is the elementary solid angle. If we decompose vector fc, on a spherical 
basis 



then 




k = Ko\ sin 61 sin (43) 



=||p||d||p||d</. (44) 



(27r)2 

^^d^^, (45) 

since fc = p + aQ{p)ez implies that ||p|| = Kq s'm9, ao(p) = -f^ocos^, and 
d^fc — sm9ddd4i. From equality 1(42(1 we find 

J^ncoh fc) = ^0 COS ^ j.ncoh (^R^p)^ (46) 

and from 1(41(1 we write the cross-section as a function of X'"'^°''(ii, fc): 

47rcos6'X™™''(i2,fc) = 5^™™''(fc|fco) : J^^'{ko) (47) 
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with KqU =p + ao{p)e^, Koko =Po+ ao(Po)^z, and 

J^Kkn) = J^^Po) , (48) 

= '""^" ^°'(Po)^^°"(Po). (49) 

We recover here the generahzation of the definition used in the scalar radiative transfer, 
where the cross-section is given by 

a"(fe|fco) = . (50) 

J{ko) 

where I{R, k) is the usual scalar specific intensity, and i7(feo) the incident power flux. 
3. Bethe-Salpeter equation and specific intensity 

To obtain the bistatic cross-section ct'"^"'*^ we are now going to calculate the specific 
intensity I by using the Green functions defined in I. We have shown, in particular, 
the following relationships: 

-^00 _ -^0 -^01 ^11 -^10 , , 

Ggy = Gg +Gs ■ Tgy ■ Gg , (52) 

where the operator T^Jv satisfies the coherent potential approximation < T^sv >v= 0. 
In using this approximation, the average tensorial products of the Green functions are 

given by 

< Gsv Gsv »sy=< Gg Gg >s 

+ <Gg(g)Gg >s:<Tg >g:<Gg (E>Gg >s, (53) 
— 11 — 11* — 11 — 11* 

< Ggy O Ggy >SV = < Gg (g) Gg >s 

ttII ^=^11 :^11* 

^<Gg®Gg >s:<Tg >g:< Gg (g)Gg >s , (54) 

where we have defined T^g =< T^gy (S> T^sv ^v- In these equations, we use the 
following convention for the tensorial product of two dyads f{r, ro) and g{r, ro): 

J (3g{r,r'\ro,r'o) =J{r,ro) (^g{r',r'o) , (55) 

and for the product between two tensors, 

[M^ :M\r,r'\ro,r'Q) = / d^rid^r[M^{r,r'\ri,r[):M'^{ri,r[\ro,r'Q). (56) 

JVi 

We now introduce the intensity operator V and the Bethe-Salpeter equation satisfied 

by < Ggy «) Ggy >y: 

— 11 — 11* — 11 — 11* 

< Ggy (g) Ggy >5y=< Gg (g) Gg >s 

+ <Gs (g) Gg >g: V :<C Ggy gi' Ggy >sy • (57) 
In the next section, we will give an approximate expression for the intensity operator 
V which will not depend on the rough surface profiles z = hi{x) and z = h2{x) but 
only on the properties of the scattcrcrs. It will appear that the tensor "P describes 
the intensity scattered by a particle taking into account the correlations with other 
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particles. To simplify the notations, we introduce two new tensors G^sv» ^^'^ Qoa<s> 
defined by 

^<<sy» ® Gly* »5v , (58) 

C<5>=<g"®g"*>s, (59) 
making the Bethe-Salpeter equation (|57|) now 

-5-11 -yll -5-11 —11 -5-11 

y ^sv» = y oc <s> + y oo<s> ■ •^<sv»- (qO) 

In definition H59|l , we have introduced the symbol oo to emphasize that the propagator 

Goo<s> between two scattering events by the particles describes either a wave 
propagating directly between the two scatterers (which is taken into account by the 

term in G^g) or a wave reflected by the boundaries (which is taken into account by 

the terms Gg in ). In iterating the Bethe-Salpeter equation (|57|l and comparing 

it with equation H54|l . we express the operator < T^g >g as a function of the intensity 
=11 

operator 7^ : 

<Ts >s =P +V ■■goo<s>-T +^ -GooKSy-T ■Qoo<s>-T +■•■. 



V +v 



pjll Tfll ttII 

^oo <S> + ^oo <S> • ^ ■yoo<S> 



^" . (61) 



The term in the bracket is identical to the right-hand side of equation (|60|l . and 
equation (|61|) is written 

< r" >5- ^" + ^" : ^<<sv» : ^" • (62) 
By introducing equation (|62|l in (|53|l . we have 

-^00 pjOO -^01 —11 -^10 -=fll —11 -p^ll —11 -^10 

(63) 

where we have defined 

G«sv» =« Gsv Gsv (64) 
^oo«sy» =< Gs ® G5 >s (65) 

7^10 -^10 -^10* 

C?<5> -< Ggy ® Gsv- >s , (66) 
TjOl —01 :^0l* 

g<s>=<Gs®Gs >s . (67) 
In equation H63() . we have expressed the intensity in the medium produced by a 
source in the medium 0, which is described by the tensor G^sv:s> ^ function of the 
tensor G<^sv^- intensity operator P represents the intensity scattered by 

one particle, equation 1)63(1 decomposes the scattering process in three terms: the first 
one where no scattering in the volume happens, the second one is the single scattering 
term, and the third term contains higher scattering contribution by particles. 

To calculate X™'^°'' in equation (|37|l . we need to connect the scattered intensity 
-^sV with the Green function Ggy. For a source of radiation produced by a current 
3 sources ^hc incident field is given by 

E°\r) =iu;iiyac d^roG^{r,ro) -jsourceiro)- (68) 

JVo 



Electromagnetic wave scattering from a random layer with rough interfaces II 10 

The field E^^y in the medium produced by the source and the waves scattered by 
the random medium and the rough surfaces leads to 

Elyir)^E^\r) + E%{r), (69) 

^iuifi^ac d^ro'^svi^^^o) ■ jsourceiro) ■ (70) 
JVo 

If we define a new function G by 

then the scattered intensity is given by 

E%'y{r)=iujn^ac d^ro'^sv°'ir,ro)-jsourceiro)- (72) 

JVo 

The tensorial product appearing in definition of the Wigner transform of the 
specific intensity 1"'-'^°^ is thus 

E°UR+^)®E%*iR^^) = (cM^ac)' 

y (27r)3 (27r)3 (-^ + 2 ' '"o) ' Jsource(^o) G-^y (-K - -, Tq) • jgQ„^^g(rQ) . 

(73) 

After some calculation and by using development l|63|l . definition 17111 . and the 
decomposition Gg = Gq + Gg^ introduced in I, in equation (|73|l . we obtain the 
three following contributions for the specific intensity defined by H37|l : 

—incoh —incoh —incoh —dncoh /^4\ 
^ —^L=0 +^L=1 +^Ladderj (74) 

with 

ITltiR,k)= eo [ d're-'"-^ [< £;°^(H + J) ® £;°^*(H - J) >s 

JVo ^ ^ 

- < El^iR+^) > ^ < El^*{R-^) >s] , (75) 



IT^f{R,k)= eo [ d3re-"=-C?<5>:P":<£;"®£;">s(il+J,i?- J), (76) 

JVo ^ ^ 

rZdLiR, fc) = eo / dVe-''^- g'^sy ■ ^" : G<<sv» 



:7'":<£;i*®£;i*>5(fi+|,H-|). (77) 



Here 



eo = ^ , (78) 

and E^'^{r), E]f{r) are, respectively, the field scattered by the boundaries without 
any interaction with the particles and the field transmitted by the boundaries inside 
the slab before any scattering by the particles: 

E%\r) = iuj ^x^ac d3roG5+""(r,ro) • j,o„^ee(^o), (79) 

JVo 

E\\r) = \u;^l,ac I dVoGs (r,ro) • j_(ro) . (80) 
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From the relationships in (|37|l and (|38|) and the decomposition 1)74(1. we can also 
decompose the incoherent scattering cross-section in three parts: 



—incoh —incoh , —incoh , —incoh 
a — (T T + (T r,-i + <T , 



' L=0 " L=l Ladder- (81) 

In section (jlOII . we will add a fourth contribution to the cross-section taking into 
account the enhanced backscattering. From equation IjTSII . we notice that the term 



^mco/i determined by the scattering matrix S 
we can write the field scattered by the rough surfaces E'^g as 

instead of equation (|79|l . The term ct™^^'' is determined by the tensor V^^ and by the 

-p;0+l± -p;l±0- . , ^ r • -p^^ 11 -1 

scattering operator o , o since the Green tunction G and the transmitted 
field E^^ appearing in the equation iTT^ are 

(r, ro) - g"/'" (r, r,) + G°/'+ (r, Vo) , (83) 



-^0+0- 



since in using the results of I, 

Os 



(82) 



-^0+ lao , , 

(^,^o) 



(27r)2 (27r)2 

0+ lao 



- — J__P_^_^0_ ip-x-ipg-aso-\-iao{p) z-aQiae{Po) zq 

x'S 



giQi(p) z ^ 



Po Po 

l + ,0 



(84) 



(85) 



where ao = ±, and K'^ — Re{Ke)- 

On the other hand, to get the term o^Ladder^ need to calculate term ^^5^-^ 
which appears in equation H77|l. As this equation is a Wigner transform, it is 

convenient to write the Bethe-Salpeter equation (j^UJl, satisfied by Q<^sv^^ 
introducing the Wigner transform of a tensor by 



g{R,k\Ro,ko) 
and we obtain 



g<<s^»(H,fe|fio,fco) =^?"<s>(-R,fc|-Ro,feo)+^ d^Hid^ila J 
xC<S>(-R: k\Ri,ki) : p"(Hi, fci|i?2, fcz) : ^<<sy»(-R2, fczlHo, fco) . 



d'^fci d^fc2 
(2^(2^ 



(87) 



In section (O, we will write explicitly the different terms of this equation under the 
ladder approximation, and in section © , we will derive the radiative transfer equation 
from it. 



4. Intensity operator V and modified ladder approximation 

In the previous section, we have introduced in a formal way the intensity operator 
V to write the Bethe-Salpeter equation. This intensity operator can be determined 
in using the energy conservation principle. In fact, by introducing the Dyson 
equation, defined in [Appendix C| and the Bethe-Salpeter equation (|57|l . in the energy 
conservation equation, we obtain a Ward identity which is a relationship between the 
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mass operator M ( [Appendix C| l and the intensity operator V . The mass operator 
is known since it is a function of the effective permittivity ee- 

M'\r,ro)^{ee-ei)S{r^ro)l. (88) 
In paper I, we have shown that under the Quasi- Cry stalhne Coherent Potential 
Approximation (QC-CPA), the scalar ee satisfies a non-linear system of equations. 
In using these equations in the Ward identity, we obtain an expression, called the 
modified ladder approximation, for the intensity operator V which satisfies the 
energy conservation ^] : 

V^\r,r'\ro,r'o)=n J d^r^ < C^V.r, (^-^o) >v-r, <c'sv,r,ir'\ro) >v-r, 

^r? j I dV, d^r, h{rj ~ n) < C^V,., (r|ro) >v-r, ^ < C^V*., (rVo) >v-r, , (89) 
where 

hir) = gir) - 1 . (90) 

The term < C^gy.^^{r\rQ) >V;rj has been defined in paper I and represents the 
transition operator for a scatterer located at rj which takes into account the correlation 
with scatterers close to the point rj through the function g{r). Hence, the first term 
in equation H89|l describes the scattering process by one scatterer located at rj, and 
the second term represents the interference process between a wave scattered at rj 
and another one scattered at r/. We do not reproduce here the derivation of the 
Ward identity and equation (|89|l since the demonstration is formally identical to the 
infinite random medium case and is well documented I34L IHl I45| . However, we 
must notice that equation H89|l is valid only for the static case when the harmonic 
dependence oj is the same in the left and right hand sides of the tensorial products of 
all equations previously written. In the dynamic case, the Ward identity contains 
a new term, taking into account the time delay that the wave undergoes during 
scattering, which modified equation (|89|l . The derivation of this dynamic Ward 
identity and the modification on the radiative transfer equation are described in 
references |2H1 HHI (47,, 48, 49, 50, 51 , 52, 53, 54. 

In I, we have introduced the following notation for the average transition operator 
of a scatterer located at the origin: 

cl\k\ko) -< CsV,..=o(fc|feo) >y;r,=0 , (91) 
where the Fourier transform of the transition operator is defined by 

<clv,r,ik\ko) >v;r,^ J d^r d^ro exp(-ifc-r+ifco-T-o) <cl^v,r,{r\ro)>V:r, -(92) 

In I, we have also shown that the transition operator for a scatterer located at point 
rj can be expressed as a function of {k\ko). 

< CsV,.,(fc|fco) >y;.,= e-'('=-'^«)-'-^ cl\k\ko) . (93) 
In introducing the Fourier transform l|92|) and the properties ()93(l in equation H89|l and 
taking the Wigner transform H8t)|) of the intensity operator V , we obtain 

V^^iR, k\Ro, ko) = nw{k — ko) 

d^K ,r^jT,_n„^^ii,K ,,K . , —11*, K ., K 



e 



(27r)3 oV2''2 ''' °'2 ' 2 

(94) 
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where w(fc — fco) is the structure factor of the medium identical to the one defined in 
scattering by X-rays [331 ■ 

wik-ko) = l + n J d-V e-'('=-'=«)'- [g(r) - 1] , (95) 

with ^(t") the pair distribution function. As the system under study is invariant along 
the X, y axis, we introduce the Fourier transform along these directions: 

v'\z,k\Zo,ko;P) = J d^(X -Xo)e-"^^-^°>^P^\R,k\Ro,ko), (96) 

with R = X + Zez, and Rq = Xq + ZqCz- From equation H94|) . we deduce that 

V^\z,k\Zo,ko;P) = nw{k - ko) 

X ldK.e^-^^'-'o)cl\^ + k\^ + ko)^cl'*{~^ + k\-^ + ko), (97) 

where Kp = P + K^Ez- Expression H97|) is rather complicated since the operator 
11 11 11 ^ 

V is non local {V ^ for Z ^ Zq), and the operators are off-shell evaluated 

since we do not have necessarily 

||ii:p/2 + fe|| = ||Kp/2 + fco||, nor || - fs:p/2 + fc|| = || - fi:p/2 + fco|| , (98) 

in equation (|97() . Accordingly, we use an on-shell approximation |30| on the transition 
-=11 

operator : 

Cl\±^ + k\±^ + ko)^ C"(fci±|fe;±) , (99) 

with 

P^ = ±y + fc± p±=±y + fcoi, (100) 

fcp± = + sgn{kz)a'^{p^)ez , fe^i ^ + sgn{koz)a'^ip^)ez , (101) 

where k± = (k-e.x)ex-'r{k-e.y)e.y, and sgTi(kz) is the sign of kz- In this way, the incident 
and scattered wave vector on the particles have the same norm | jfcpi 1 1 = | |A;^± 1 1 ~ K'^. 

With this on-shell approximation, the intensity operator V becomes localized at the 
point Z = Zq, and we have 

p'\z,k\Zo,ko;P) = 5{Z-Zo)v'\k\ko;P), (102) 
p'\k\k„;P) = nw{k-ko)cl\kl,+ \kl^)<E>cl^*{kl-\kl-). (103) 

We must also emphasize that the correlations between the scatterers appear not only 
in the structure factor w{k~ko) but also in the scattering operator , which satisfies 

/i3 T 
h{k - k,)tl\k\ki) ■ Gr(fci) • c"(fei|fco) , (104) 

with 

h{k - fci) = J d^r e-'('=-'=i) '- [g{r) - 1] . (105) 

Here, t^^ is the scattering operator for a particle located at the origin which, contrary 
to Cg , does not take into account the correlation with the other particles. 
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5. Bethe-Salpeter equation and ladder approximation for the rough 
surfaces 

As discussed in section we need to know the tensor 

^<<sy» (g) Cfy :^sv (106) 

to determine the specific intensity X™^'^^'^^^, according to equation l|77|l . We know that 
this tensor verifies the integral Bethe-Salpeter equation H87f) : 

/f d"^fc d'^Ai 
^ d^Hid^ila J ^^^^3 ^^^^l 

xgl^^s>iR^k\Ri,ki) ■.V^\Ri,ki\R2,k2) : g<<sy»(i?2, fe2|fio, ^o) . (107) 

Between two scattering processes on the rough surfaces, the wave interacts with the 

particles, and accordingly, it is reasonable to admit that the scattering events on the 

— 11 — 11 — 11* 

rough surfaces are uncorrelated. Hence, the propagator G^sy ^ >s in 

equation (|107|l is approximated by 

< Gg (g) G5 >s= Gi (g) Gi + < G5 (g G5 >s + <Gs g) G5 >s 
+ < g]^^~ (g gI^^^ * >s + < gI^^^ (g gI^^^ * >s (108) 
since the Green function G5 has a development in five components: 

Gl' = gT + Gr^ + Gr^ + G^- + Gl-'- . (109) 
In using condensed notation, equation (|108|) is written 

C<5>=C+ E sTst\ (110) 

a,ao — i 

g'J^^gT^gT*, (111) 

gTs7=<GT''"^GT"*>s,, (112) 

^"<s>=<g"®g"*>5, (113) 

and a, oq are the sign + or — . Furthermore, the Green functions Gg " are defined 
by using the scattering operator b . iiach operator b can be decomposed 
by using the operators R and R describing, respectively, the scattering by the 
upper and the lower rough surface. Then, the hypothesis of non-correlated diffusion 
on the rough surfaces must also be applied on each term of the operators S " °° . For 
example, the Green function Gg depends on S , and this operator has the 
following development: 

S'-'' = R°' ■r"''+^'-r"''-r''-r"'' + .... (114) 

i_i_ i_i_ 1—1—* 

With our hypothesis, the tensorial product ^^5^ ^<Gg g) Gg > 5 contains 

the following terms: 

—1-1- — I-I-* —01 T^Ol* -FTff21 —H21* 

<S (gS >s=< -R <g-R >s-<R iXi-R >s 

-=pi TiirOl* -^H21 —H21* —01 -=01* -p^H21 —H 21 * 

+ <R g)i2 >s-<R g)-R >s-<R <^ R >s-<R g) -R >s 
+ .... (115) 



with 
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Under these approximations, that corresponds to the ladder approximation for the 
rough surface contributions; the Bethe-Salpeter equation H1U7|) describes the scattering 
processes depicted in Figure ^ where the waves on the left and right sides of the 
tensorial products follow the same path. As shown in Figure 13 we can also use 




Figure 1. Scattering processes under tlie ladder approximation. The solid line 
and the dashed line represent, respectively, the wave on the left and right hand 
side of tensorial products such as Q ^sv^ =<KGsV®''sV*2*sv- 

Feynman diagrams pi I24L 1^ to represent the different terms q'^Q^sv^ the ladder 
approximation obtained by iteration of the Bethe-Salpeter equation 1)107(1 : 

+ C<5> :C<5> :^"<S> + --- • (US) 

We now examine in detail the different terms in the Bethe-Salpeter equation (|107() . 
Since, by hypothesis, the medium and the surfaces are statistically homogeneous, we 
introduce a Fourier transform along the X, Y axis: 

G{Z, k\Zo, ko; P) - / d'(^ - ^o) e-'^^-^^' -P g{R, k\Ro, ko) , (117) 
with 

R = Xe^ + Yey + Ze:, = X + Ze^ , (118) 
Ro = Xoe^ + YoBy + Zoe, = Xq + Zoe, , (119) 

and where y is either y<gsv»' »oo <s>' »<s> i »oo- With this Fourier transform, 
equation (|110|l is now 

g]^^s>{Z,k\Zo,ko;P)^g]^{Z,k\Zo,ko;P)+ G<s>'iZ,k\Zo,ko;P) . (120) 

a,ao — i 
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y«s 



Gi 



?11 

'oo <S> 



Gi 



Gc 



Gi 



Gi 



Figure 2. Diagrammatic representation of tlie ladder approximation. 



Furthermore, in calculating the terms G^s> ^-nd in equation (|120|l and using 
the Weyl representation of and the Green function decomposition in terms of 
the scattering operator g^""^""" given in I and mentioned at the end of section (PJ, we 
obtain the following decomposition: 

-lalao, . . , ^2cA ae(P+) +ae(P~) 

y<S> {^^'^\^o,l«o;P) = [2n) dlk^-a 

x<5 Uo.-ao^^^^^^^— g^s> {Z,p\Zo,Po;P), (121) 



a,ao — ± 



ae{p+) + alip- 



xd[ko^-ao [Z,p\Zo,Po;P) , 



(122) 



(123) 



where a and ao are the signs + or — , k = p + e^, fco = Po + ^oz ^2, and 

+ P + P 

P =P^Y' Po=Po±y- 

Under the quasi-uniform field approximation ^21) we approximate ae(Po') + a*(p|^) 
by ae{p) + a*(p) = 2Reae(p), and equations p21l) and p22|) become 

g<i>''(Z,fc|Zo,feo;P) 



(27r)2^(fc2 - aa^(p))J(fco. - ao a'M)g'^s> iZ,p\Zo,P,; P) , (124) 



^11 



g^{Z,k\Zo,ko;P) 

= X! (27r)^'5(fc2 - aae(p))'5(fco2 - aoae(Po))^'i^^"''(^>Pl^o,Po;-P) ^ 

a,ao— i 
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(125) 

where a'^{p) = Re[ae(p)]. The Dhac distributions insure that the wave vector 
directions k and fep have constant norms ||A;|| ~ K'^, ||fco|| ~ K'^. Furthermore, 
we found the foUowing expression for the tensors G<^s> ^^"^ ^oo 

4ae(pJ)a*(Po) 

Gc^ ''\Z,p\Zo,Po;P) = (27r)2(5(p-po)7i°(p+)(g)7i°*(p")(5a,ao^a.ssn(Z-Zo) 

giaAQ, (p:P) {Z-Z„) 



4ae(p+)aj(p-) 

where 



(127) 



r _/l if a = sgn(Z - Zo) , . 

da,sgn(Z-Zo) - | g if a ^ Sgn(Z - Zq) , ^^^^^ 

with sgn(Z — Zq) the sign oi Z — Zq, and 

Aa,{p;P) = a,{p+)-a:{p-), (129) 
7i"(p) = (7 - k;k';) = + eH{p)eH{p) , (130) 

kl^ =p + aa,{p)e,, fcp" = fc^Vllfci'^ll . (131) 

The tensors 5^ describe the intensity scattered by the boundaries and are defined 

by 

{'2n?5{0)t;':s%{p\p,-P,) ^KS'""'"" ^s"''^"*{p-\Po) >s . (132) 

In introducing properties (|124|l and (|125|l in the Bethe-Salpeter equation p(J7|l . we 
demonstrate that G<^sv^ verifies a development similar to equations (|124|l and (|125|) : 

Gl^sv»iZ,p\Zo,Pa;P) 

= J2 (2^)2<5(fc,-aa;(p))<5(fco.-aoa'e(Po))^<<sy»(^-Pl^o,Po;^)- (133) 

a.ao — i 

Accordingly, we can write equation (ll(J7|) for each sign a and ao, and we have 

g:,:^;V(z,fc|z„,fc„;P) = aL"5>(^,pi^o,P„;p)+E /#ii^r 

aTt^^ (2^)M2^)^7_^ 

xaoo<5>(^,Pl^2i,P2;P):P \p^\p,;P):g^^s;^yiZ2i,p,\Zo,Po;P), (134) 
with 

^L"<5V(^,Pl^o,Po;i") =^L"'""(^,Pl^o,Po;J') +^<s> (135) 

— Ia2lai , . /rl<i2 f I'll \ -^11 /,- I lai , -=11 * , r la2 , f lai . 

P (p2|Pi;P)=riu;(fcp^ -fepjC„ (fcp+|fcp+)®C„ (fcp-|fcp^) (136) 
p± - ± y , = P2 ± y , (137) 

fep" = p + aa^(p) . (138) 

The numerical solution of integral equation (|134|l is a difficult task, and in the next 
section, we will give a differential form of this equation more appropriate in this case. 
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However, equation (|134|) is well suited if we want to obtain an iterative solution for 
^<KSv'>- example, in introducing the first term of this iteration, 

gl:];^^iz,p\Zo,Po;P)^GT"iz,p\Zo,Po-,P)+G':s>i^ 

in equation l|77|) . we obtain the double scattering contribution by the particles. The 
term Qi describes the direct propagation of the wave between the two scatterers, 
and G ^s> describes the interaction with the boundaries between the two scattering 
processes by the particles. 



6. Radiative transfer equation and boundaries conditions 

To derive the radiative transfer equation verified by ^^5y^(Z,p|Zo,Po; -P); must 
first notice the following properties of Goo ° ^^^^ G^s> ' 

+ ia Aae(p; P)gI^'"°{Z,p\Zo,Po; P) (139) 
^g<s> (Z,p|Zo,Po;-P) =iaAa,(p;P)^<s> (140) 
with 

tl^{p:P)=tl{p+)®tr{p-). (141) 
p±=p±y. (142) 

To obtain equation p39|) . we have to introduce in equation p27|) the identity 

'5a,sg„(z-z„) = '5a,+e(Z - Zo) + 5a^-Q{Z^ ~ Z) , (143) 

with O as the Heaviside function, and then use the property: 

^(z) = 5(z). (144) 
az 

Hence, in taking the derivative of the Bethe-Salpeter equation H134() . we obtain the 
radiative transfer equation satisfied by ^^5.y~^(Z,p|Zo,Po; P) H- 

,7^,7 ^ rj. (27r)^ (5(p - Pp) -la 

^G^^Sv>>{ZMZ..P,.P)^j^^j-^^^^I^^{p,P)5a,a.a5{Z-Z,) 



Aae{p+)al{p , ^^^^ 

+ i a Aae (P ; P) ^<<5y » {Z ,p\Zo,p^; P) . (146) 
To solve this differential equation, we need boundary conditions on the interface. In I, 
we have shown that the scattering operator S can be decomposed as a function 

1 To get this equation, we have supposed that the operator c]^{k\ko), which defined the operator 
■p^^ , is transverse to the propagation direction defined by fc: fc ■ c]^ (fc|fco) = 0. In this case, we have: 

This hypothesis is satisfied in the far-field approximation. 
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of the scattering operators il*^^ and R^^^ of the upper and lower rough surfaces. In 
using these equations, we easily show the following relationships between the scattering 
operators: 





-1- 




(147) 


s'- 


-1+ 




(148) 


s'- 


fl- 




(149) 


s'- 


fl+ 


^r"''-s'-'\ 


(150) 



where we have used the notation, 

[f-g]{p\p,)^ J^f(p\p^).g{p^\p^). (151) 

With equations H126I127I147I150|) and the definition of R given in I: 
we obtain 

^L"<s>(^ = 0,p|Zo,p„;P) 



/ §^T^°^<s>iP\Pi-.P)--G]^<s>iZ^O,P^\Zo,Po;P), (153) 
^Zo,Po;P) 

^ n^^^s>iP\Pi:P) ■■ ^L"<s>(^ = -H,p,\Zo,Po; P) , (154) 



G'^Xiz = -H,p\Zo,Po;P 

12 

{2n) 

where we have defined 

{27r)H{0)'R°^^s>iPW.P) =<lf\p+\p^)®^\p-\Po) >S, (155) 

{2n)H{0)nl\s>{p\Po;P) ^<R'\p+\pr)^R'\p-\Po ) >S ■ (156) 

Then, with the help of the Bethe-Salpeter equation (|134|l . we easily demonstrate that 

Q ^gy~^{Z,p\ZQ,pQ; P) verifies the following boundary conditions on the upper and 
lower rough surfaces: 

G<<sv»iZ = 0,p|Zo,Po; P) (157) 
= / ^K<s>iP\Pi-, P) ■■ G\<s7»iZ = 0,p,\Zo,Po; P) , (158) 

K<sv»iZ = -H,p\Zo,Po; P) (159) 
- / ^Kj<S>iP\Pi;P) ■■G^KSv'yyiZ ^ -H,p,\Zo,Po;P) ,{m 

With equation 114(j|) and the boundary conditions l|158|) , l|16U|) we thus have obtained 

a closed system of equations sufHcient to calculate the tensors G]^^sv»- section (|12() . 
we will show how to rewrite the radiative transfer equation l|146f) in a more classical 
form. 
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7. Calculation of wf^f" 

In sections Q and ((SJ, we have shown that the first term of the incoherent scattering 
cross-section is given by 

^^Xrro^i?, k) = i2n) 5ik^ - ao(p)) ^n\p\Po) ■ J"\Po) , (161) 

TT 

with k = p + k^e^ and 

- < E"s'{R+'^) > ® < ^) >5] . (162) 

Furthermore, we know that the scattering intensity E^'^ can be expressed as a function 
of the scattermg operator b by 

E°six,z)^ J (^e'f-+'""(f)^5°+'"-(p|po)-i^°''(Po)- (163) 
Then, we find 

crL"=o (pIPo) = ^^«5®<s> (p|Po)> (164) 

TT 

, -jiO+O-mco/i . J „ , , 

where t>^^gy is defined by 

{2.f5{0)^:rsr\p\P.) = < s"^°'{p\Po)^s"^'-{p\Po) >s 

- < S"^"- (pIpo) >s ® < S"^"' (pIpo) >s . (165) 

8. Calculation of ff^^if 

For the first order scattering contribution by the particles, the scattered specific 
intensity is 

Xt=t{R.k) = eo / d^re-''" '- < ® G^'* >s: :< -E^* ® i^^** >s , (166) 
and the cross-section verifies the fohowing equation: 

^^r£Zt{R.k) = {2n)S{k^ - ao(p))ar=ff (pIPo) : J"'iPo) ■ (167) 

TT 

In introducing equations (|83l 184(1 and IHHwhich connect the Green function and 

the transmitted field E^* to the scattering operators S and S , in ((166(1 we 
obtain 

— o.. , ^ ^ Mp)! V /■ 1 (e'[^^A..(p,)-a,Aa.(p,)]H_i) 

x5°^''<s>(p|P2; - 0) : P''^^'''^ P = 0) : ^"^^^(pilpo; P = 0) , (168) 

where 



{2n 


2^(0)5°+<';;(p|Po;Po) 


=<s°+"^^ (p+ 




s (p 




(169) 


(27r 


2<5(0)4'^^^°;(p|Po;Po) 


=<s"'''-{p+ 




S (p 




(170) 




^'"^(P2iPi) = nz«(fc;,;^- 


K:)cl\k';^ 


lfe'?) 

Pi 


^G"(fc^"_= 
" P2 




(171) 
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with 

+ P + P + P + P , , 

P -P±Y' = P2 ± y . Pf = Pi ± y , Pq ^Po±Y' 

kl'' =p + aa',{p)e,, (173) 

and 

Aae(p)=ae(p)-K(pr (174) 

= 2ilm[ae(p)]. (175) 



9. Calculation of 



Similarly to the previous section, the multiple scattering contribution is given by 
T£aMer{R. k) ^ c, f dV e-'^- < cf * >s: : ^<<5y» 

JVo 

■.V^^ :< ^Sy* >s , (176) 
and the cross-section is 

^IZ^^^AR, fe) = (27r) <5(fc, - ao(p))5^ra.l.(p|Po) : J'\Po) ■ (177) 

TT 

From equations JHSH^EHl) and properties (|117lll33|l of the tensor ^^sv'»' ^^'^'^ 

-=vncoh ^(P)! \- C A7 A7 [ ^'Pl ^'P2 d^Pg d?p^ 

f^LaddeAPlPo) - }. / dZ43ClZ21 / 



(27r)2 (27r)2 {2nf (2tt)^ 

1(3-4 /I rt ;5^1'^4lfi3 



ai ,a2 ,03 ,a4 — i 

-ia4 Aae(p^) Z43+iai Aae(Pi) ^21 

p 5;7s;(p|p4; P = 0) : 7'^"^^"^(P4|P3; = 0) 



4|ae(P4)l' 

:^43,P3l^2i,P2;J' = 0) : r'''''''\p,\p,:P = 0) : 5;"^°;(pi|po; P - 0) . 

(178) 



10. Enhanced backscattering and reciprocity 

During the last decade, several studies have been concerned with the enhanced 
backscattering Ii[T3ElE3|2niE7||nniEIlE3EHlEiEniinilE^ This phenomenon 
produces a peak in the backscattering direction (fc = — fco) due to the interference 
of waves following the same path but in opposite directions (Figure (|2Jl)- The most 
elegant way to introduce the enhanced backscattering in the theory is to use the 
reciprocity principle |63[ I64| . Heuristically, it means that when we exchange the 
source and the detector position (respectively given by Vq and r) and also their 
polarization, the field measured is the same. In using our Green functions, we translate 
this statement by 

G'svir,ro)^[G"sy{ro,r)f, (179) 
where the transpose of dyad /, which exchanges the polarization, is defined by 

— T — 00 

/ ■ E ^ E ■ f. If we decompose the dyad Ggy{r, Vq) in the fixed basis [e^, ey, e^], 
GZ{r,ro)^ J2 Gfy{r,rohe,e,, (180) 
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Figure 3. Scattering processes which are taken into account by the most-crossed 
contribution. 



we re-write equation (|179ll under the following form: 

Gfy{r,roh^Gfy{ro,r),,. (181) 

The properties in (| 179(1 can be demonstrated by using the Green theorem |65| . By 
applying these properties on each element of the tensor, 

75-00 , , . , . -^00 , , -p=00 * , , , , 

^«5y»('^. r ko, T-o) Gsvir, ro) (E) Ggy {r , r^) >5y , 



-7=^00 -^00 * , ,1 , , r-=00 -^00 * , , I 

Gsv » Gsv (r, r |ro, r^) = [Ggv » Ggy (ro: r^r, r 
«) Gsv (r, r \ro, ro) = [Ggv ® Ggy (r, rglro, r 



'MTr, 



=00 



=00* , 



we obtain three identities that must be satisfied: 

(182) 

rgy ^ ^j^gy yr , r |v-o, '-o; = I'^sv ^ ^sv v''r''or''o, '- ;j - , (183) 

^00 -^00 * / / I /X r"^00 -;=^00 /, I \->Tr A\ 

Gsv ® Gsy (r, r'Iro, = [Ggy ® Ggy (ro, r'|r, r'^)f^ , (184) 

where we have defined three transposes of a tensor Ai with the following decomposition 
(see [Appendix A| ) : 

M ^ ^ Mri'-jj' (e^ (g) Bi') {e-j (g) e^-) , 

by 

^^'^ ^Mjf.,u' , (185) 

- (186) 

(187) 

We can easily check that if two of the conditions p82ll84|) are satisfied, the third 
condition comes automatically. These properties can be reformulated by using the 



Electromagnetic wave scattering from a random layer with rough interfaces II 23 

tensor r" =< T^V ® T'sV >v- 

Tsir,r'\ro,r',) = [r"(ro, r[,|r, r')]^- , (188) 

r"(r,r'|ro,r;,) = [Ts ir,r',\ro.r')f- , (189) 

r"(r,r'|ro,r-;,) = [r"(ro, r'|r, rj,)]^- , (190) 

In tact, it we suppose that the scattering operators .5 , b , and o verity the 
reciprocity principle, it can be easily checked that the Green functions G , G , and 
G verify a reciprocity condition similar to property p79|) . Now, in section ||2Jl we 
have demonstrated that 

Gsv Gsv =Gs (^Gs +Gs (E)Gs : : , (191) 

and from the reciprocity of the Green functions g"", G ^ , G^ , we deduce properties 
H188ll9Uf) . The first condition (|188|l is a reciprocity condition on energy. If we write 
it using a Wigner transform, we have 

Tl\R,k\Ro,ko) = [Tl\Ro,^ko\R,-k)f^'' . (192) 

Equation (|192(l signifies that if the source and detector positions are exchanged 
(R <-!■ Rq) with their polarization (transpose T^r) and if the incident and scattered 
wave directions are exchanged and inversed {k ^ — fep), then intensity measured is 

the same. In section 0, we have also demonstrated that the tensor Tg is given by 
equation (|^ : 

rs{R,k\Ro,ko)^P iR,k\Ro,ko)+ d'Rid'R2 —^—^V (i?, fe|i?2, fea) 

: <Gly®Gl^y >v {R2,k2\Ri,ki) ■.V^\Ri,ki\Ro,ko). (193) 

Furthermore, with the help of the Wigner representation (|94ll of V^^ and the Bethe- 
Salpeter equation (|107|) . we demonstrate that 

v'\R,k\Ro,ko) = [p"(Ho,-fco|-R,-fc)]^"«, (194) 
<Gly(g)G'sy >y {R,k\Ro,ko) = [<Gly(g)G''sy >y iRo,-ko\R,-k)f^'' , (195) 

if we suppose that the Green function Gg and the scattering operator {k\ko) verify 
the reciprocity principle. For the operator G^ , this principle can be formulated as 

cl\k\ko)^[cl\-ko\-k)f. (196) 

In I, we have shown that this property is effectively satisfied since the scattering 
operator for one particle verifies a similar property [SHI IHZl EH ■ Accordingly, from 
properties H194|l and (|195(l and equation (|193|l . we demonstrate that Tg {R, k\Ro, kg) 
satisfies condition H192(l . 

However, tensor Tg (i?, fc|i?o, ^o) does not satisfy conditions H189|) and H19U|) . In 
fact, these conditions describe the scattering process where two waves follow the same 
path but in opposite directions as shown in figure 0. If we decompose tensor Tg 
under the form: 

T'g =r'' +Tl\^,,,,, (197) 



with 



Ladder :< G gy G gy >v: V , (198) 
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it can be checked that only T^^^^adrfer does not verify equations H189|l and (|190|) . 
In fact, we show, in using the inverse Wigner transform V (r,r'|ro,ro) of the 
tensor V {R,k\RQ,ko), that all the conditions f (|188|l - H190() ') are satisfied by 
V^\r,r'\ro,r'o). 

When we applied the reciprocity condition (|189|l to the tensor Ladden obtain 
a new tensor crossed- 

^s'crossed{r,r'\ro,r'o) = [Tl\adder{r,r'a\ro,r')f'' . (199) 
To insure that Tg verifies all the reciprocity conditions, we add the contribution 

Crossed tO Tg I 

+ Ladder + '^S Crossed ■ (200) 

This new contribution can be represented with the help of diagrams as the set of 
"most-crossed diagrams" ^H^ESHSU], as shown on Figure 0] In using the identity 



Figure 4. "Most-crossed"' diagrams. 

p99|) in the calculation of the ladder approximation, we obtain the contribution of 

S Crossed 



-ijicoh 

Crosse^ 



MP,) = — - — _____ E__^^J_„ dz„dz„y ____ 



iae{^ + P4) a*{-^ + P4) 

-las 1 02 - \ r7 T->\ ^^'-'■^ Ic^l / I 7^1 '^i 0- 



Tr 



■.G<<SV»(.Zi3,Ps\Z2i,P2;P):V - (p2|Pi;P):5^<s>(Pi|-Ap;P)J , (201) 

where we have used notations p69ll71f) and defined: 

P=P+Po. (202) 

Ap = (203) 

Aa{p; P) = a,{p + P/2) - a^p ^ P/2) . (204) 

The definition of the right transpose '^'^ when the tensors are expressed in the basis 
[e^ (p), enip)] and not in the the fixed basis [e.x,ey,ez] is defined in [Appendix D| 
It is worth mentioning that to calculate the contribution, we need to know the tensor 
in the ladder approximation, and thus, we can use the same radiative 
transfer equation to calculate the ladder and the crossed terms. The difference 
between the two calculations is the value of P which is the null vector for the ladder 
approximation and is given by equation H202|) for the crossed contribution. 
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,,aiZ,,p,;P)^iaiAa,{p^;P)l'{Zi,p,;P). (206) 



11. Another formulation 

In the previous development, the main quantity that has to be determined in 

using the radiative transfer equation is the tensor Q ^gy^{Z43,p^\Z2i,P2] P)- The 
main advantage of this formulation is that the ladder contribution and the crossed 

— 1 0,3 1(2.2 

contribution can be evaluated as soon as Q ^gy^{Z43,p^\Z2i,P2', P) is known. 
Nevertheless, this approach is rather unusual, and in this section, we introduce the 
specific intensity X(i2, p), and we show that this specific intensity satisfies the classical 
radiative transfer equation. First, we introduce the reduced intensity by 

T];UZi,P^;P) = K:<sliPi\ - Ap; P) e"^^ ^-(p^^^)^^ . (205) 

It can be easily shown that this intensity satisfies the radiative transfer equation H14t)|l 
without the source and the scattering terms: 

d —lai 

dzi 

We now define the tensor j]^{Z,p) by 

J^{Z,p)= J2 / 7:r%^ ^(p|pi;P):X^;,,(Z,Pi;P) 

y (2#(2$(2^^ (p|P3;P) : W>(^,P3l^2i,P2;P) 

:P'"^'°^(P2|Pi;i') :X;";,,(Z2i,Pi;P), (207) 
and the diffuse specific intensity X^^(Z,p) by 

t;^[z,p)=Y^ r dzi / ^^L"<5>(^,pi^i,Pi;^)::7«'(^i,pi), (208) 

where I^^{Z,p) corresponds to upward wave propagating in the direction k = 
p + a'g{p)ez, and X^^(Z,p) to down-going wave propagating in the direction k — 

p — Q;g(p)e^. We notice that in the definition of J'^ {Z,p), we have added the first 
order scattering term to the high-order scattering contribution. With these definitions, 
we can express the scattering cross-section as a function of I^^{Z,p). First, in 
noticing the following properties of the scattering operators S^^^^: 

^0+1- ^ -01 -1+1- -0+1+ ^ ^ . ^1+1+ (209) 

that we rewrite in one equation: 



s (pIpo) = j (pI^'i) 



{27rrS{p,-po)S,,,+ l{^ip,) + s'^""'{pM)] ■ (210) 

In using our hypothesis of statistical independence of scattering processes on the 
boundaries, we have 



S^<S>iP\PQ-^P) J J^^^<S>iP\Pl-^P) 



ii2nf6{p,-p,)6ao,+ t^{Po)^I±\po)+KTs>iPi\Po-,P)] , (211) 
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where 

(27r)25(0)r^<s>(p|Po;P) =<T"\p+\p+)^¥'*{p-\p^) >s, (212) 

with = p ± P/2 and p^ = pg ± P/2. Then, from definition H135() ofQ^^^gy and 
equations H126() and H127I) . we obtain 

= /S55l>'<s>(PlPi;-P):5^«>(2 = 0,Pi|Z4.P4;P). (213) 

with p^ — p^± P /2. In comparing equations pt)8l 11781 1201|l with definitions (|205l 
EnminHi) and expression (|213|l . we obtain the foUowing expression for the scattering 
cross-sections: 

-=incoh(^\^ \ ,—incoh f^\„ ^ '^ojp) 

f^L = l (PlPo) + LaddeAPWa) - 

TT 

(2^^®<S>(P|P4; = 0) :X^+,(Z - 0,P4; P - 0) , 

(214) 

—incoh \ ao(p)ao(Po) 
f^CrossediPlPQ) = Z 



/ §^'^°^<s>(^p\P4-^P) --^miz = o,P4;-P)]^" - ^^5^r=f (p|Po;-P)(2i5) 

with P = p + and Ap = (p - Po)/2. The term o^™''(p|pq; P) is obtained 
from o'^"^i'^(p|Po) in replacing in equation H168|) . P = ObyP = p+Po and Aae{p) 
by AQ;e(p;P) defined in H2U4|I . To get the radiative transfer equation satisfied by 

I °^{Z,p; P), we first derive equation (|208|) by using H139II140|I : 



^{Z,p; P) = ^^-^jl^iZ,p; P) + iaAae(p; P)X^°,(Z,p; P) . (216) 

Next, we express the tensor J'^ as a function of X^^. From definition H207|l and the 
integral equation l|134|l satisfied by G<^sv^ ; obtain an integral equation on J"^ : 

r:{z,p)= ^v''''''\p\Pr;p):t:uz,p,;p) 

a,X^±J~H J (2^) (2^) 



■.j];\Z,,p,;P), (217) 
and introducing definition (|208|) in the second term of right hand side of equation 
H217|l , we express J'^ as a function of I^^ : 



T;{Z,p) = ^v'^"'\p\p,;P):t:UZ.P^;P) 

T.J ■■^miZ,P^.P)- (218) 



ai— ± ' 
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By introducing equation ()218|l in H21fcifl . we obtain a radiative transfer equation on 



-^(g)d- 
d_^la 

dz 



-=la , _ , Oa V — ^ f d^T3i :=^lalai , , -=la-\ , _ , 



+iaAae(p; P)t'^AZ,P; P) + sl^{Z,p- P) , (219) 

—la . 

where is the source term: 

ttIi , 0-4 X — ^ f d"^/?! -=-lalai , , „, — lai , _ , , 

("'^'^^ = a,{p+)at{p-) (4^^ ^^'^^'^^ :X«„,(Z,Pi;P) (220) 

From definition l|2()8|l and the properties H1531 1154(1 , we easily obtain the boundary 
conditions satisfied byX^^j: 

^^'<5>(p|Pi;P) :X^+(Z = 0,Pi;P), (221) 

^ <<5> (pIPi ; P) : X;",(^ = -H,p,:P), (222) 

We also notice that from equation (|218|l and definition H208() . we obtain the integral 
formulation of the radiative transfer equation 1219() and it's boundary conditions 

(|22ll), 1223): 

=ltt / ^ X — ^ f , ^ f d^Pi d'^Po — lalao , „ , _ ^-Idolfli / , 

I^,iZ,p;P)^Y. dZ2i -^-^g^JyiZ,p\Z2i,p^;P):V \P2\P1; P) 

■.t:UZ2uP^;P)+ E f dZ2J^^G'^':s,iZ,p\Z2^,p,;P) 

■.V"''"'\P2W,P) ■.TlZiZ2i,P^;P). (223) 

In the next section, we will show that this formulation is identical to the one 
usually used in the phenomenological radiative transfer theory. 



12. Simplifications 

So far, we have used the two-dimensional vectors p to described the direction of 
propagation. To recover the usual formulation of the radiative transfer theory, we 
have to introduced the three-dimensional normalized vector k. For a vector in the 
medium 1 such that ||p|| < K'^ (hence, we neglect the evanescent waves), we have the 
following relationship between p and k: 

K'^k^p±a'^{p)e, or p = K'^k^ with ±a'^{p) ^ K'/k, , (224) 

where the -|- and — signs correspond, respectively, to upward and downward waves, 
and k± = kxe^ + kyBy is the projection of the vector k on the plane (Oxy). To 
achieve the transformation of p into fc, we have, in particular, to express all the factors 
involving ae(p ± y)j occurring in the previous equations, as a function of a'^{p)- We 
remark that P is different from zero only for the enhanced backscattering contribution. 
Furthermore, this contribution is important due to the exponential factor in equation 
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H2U1|) . only closed to the backscattering direction where P ~ 0. Thus, in a good 

p 

2 



approximation, we can use a Taylor development of ofg (p ± ^ ) , and we have 



ae(p±y)~a,(p)±y • Vae(p), (225) 
= ae(p)±;5^, (226) 



since ae{p) = — . Moreover, in most cases, we have K'J <C K'^ for 

Ke — ifg + ^K'J, and a Taylor development of ae{p) which gives us 

«e(p)««l(p)+i:^, (227) 

with -\/(ReKep"^^p2 ^ RcY/i^e — = «e(p)- particular, from equations H226|l 
and H227(l we have 

. . X K'K'J P-p , , 

In using these approximations in the radiative transfer equation (|146|1 . we obtain 
— lip-P + oaJp)— 1 g^sv:s>{Z,p\Zo,Po;P) + Keg^sv:^{Z,p\Zo,Po) 

(27r)^ (5(p - Pq) -la, , -=la, , t.,^ ^ , 

= iK^a'Ap) ^ ^^^^^ (p)<5a,a„a<5(Z-Zo) 

1 V -V /" d^Pj^ — lalai , . , — laiao 



4i^>^(p) ^^^y (27r) 



ai— ± 

with — 2A^g , and 



E /^^'""(piPi)^^"r;»(^,Pii^o,Po;p), (229) 



-=^lalai , , , -=^lalai , , „ , , , 

V \p\p,)^r \p\p,;P^O). (230) 

In writing P = in the tensors V ^ (pIPi! P) and /|'|_'^(p + P/2)§5/|'|_"(p — P/2), we 
impose that in the enhanced backscattering contribution, the direction of propagation 
is the same for the two waves which we consider to calculate the intensity. Hence, 
the contribution of P appears only in the first term of equation H229I) . In the integral 
form of the radiative transfer equation given by equation (|134|) . these simplifications 

allow us to write the Green functions G]^s'> ^^'^ ^m^"" given by p2t)l 112711 as 

—lala„ (27r)^ 5(p - Pq)-1q T^la*/ ^r. r 

Goo (Z,P\Zq,Po;P) = —. — — (Po) ® -fj. {Po)^a,ao6a,sgniZ-Zo) 

4 \ ae{PQ)\ 

pylalao 5g<5>(p|Po;P = 0) 

y<s> [Z,p\Zo,Pq;P) = - 



2 



4|ae(Po)l 

X exp (a (-Ke + i^T^ ) Z - ao (-Ke + i^^^r^ ) ) , (232) 

where the factor P appears only in the exponential terms. 
Furthermore, if we remember that 

a:,:^;V(H,p|iio,Po) - / 0e'^ (^-^")^:,:^;V(Z,p|Zo,Po;P), (233) 
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with R = X + Ze^, and Rq = Xq + ZqE^, then equation H229|l is now 



(2^)2<5(p-p„)_i. 



4i^X(p) 

1 



-f± (p) «> -f± (p) (5 S{R~Ro) 



E 

ai— ± 



(27r)2 ^ 



—lai Ian 

«;sy> 



(-R,Pl|-Ro,Po) : 



(234) 



with kp"" — p + aa'^{p)e.z and fcp = k^p /K'^. 

There is also a another way to describe a vector k of constant norm if^ other 
than using the decomposition fcp" = p + a ae(p) ^i^h a = sgnikz)- We can introduce 




Figure 5. Wave vector decompositions. 9 G [0, tt], and G [0, 27r]. 

the normahzed vector k such that k = ||fe|| fe, and 
/ sin cos 6* \ 

k = 




(235) 



In this way, for every function F{k) which contains a Dirac distribution fixing the norm 
of the vector fc, we can define two new function F{k) and F°-{p) (that we differentiate 
by their variable k or p) by 



F{k)=^^5{\\k\\~K)F{k), 



such that 



F{k) = - a«e(p)) F^iP) , 

a=± 

f ^ Fik) = / d'k Fik) 7^^ ^"(f 



(236) 
(237) 

) , (238) 



where d^fc = sin 



is the integration on the solid angle. We can connect the 



functions F{k) and F{p) from equation (|238|l by 

d^p 



F''{k)d^k = F^ip) 



(27r)2 ' 



(239) 
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where we have defined for a = ± 

F''{k) = Q{ak,)F{k). (240) 

with the Heaviside function. From Figure [S] we have 

\\p\\^Ks\ne, (241) 
a',{p) = K\cose\, (242) 

since for a wave in the medium 1, we have ||fe|| ~ K'^. Then 
d^p 1 



(27r)2 (27r)2 



p||d||p||d0, (243) 
dfc , (244) 



F''{k) = ^ F'^(p) . (245) 



(27r)2 

so that equation (|239|l can be transformed into 

{2n) 

In particular, for the Dirac distribution F{k) = (27r)^ S{k — ko), we have 

F{k) ^ b(k - feo) , (246) 
^^(p) = (27r)2<5(p-p„), (247) 
and from equation (|245|l 

J(fc-fco) -if^«^(p)<5(p-Po)- (248) 

Hence, rather than using the variable p and Pq for ^<sv'>, we can introduce the 
vectors k and fco and, defined similarly to equation H239() . new tensors such that 

y«sv»(-R.fe|-Ro,Ko) = (27r)2 ^«g^»(-^'Pl-^0'Po) j^27r)2 ' 

with A; = K'^k = p + aa'^{p)e.z and fco = K'^k^ = Pg + '^oQ!e(Po)ez. With this 
definition, the radiative transfer equation H234() is now 

• vy^sy3,(if,fe|i£o,fco) + K^g ^sy^{R,k\Ro,kQ) 

= \, Sa.ao SiR ~ Ro) I± (fc) /_!_ (fe) 

+ Y. I d2fci^'"'"'(fc|fci):^<<5y>>(^,fei|^o,fco), (250) 
where a = sgn(fc • e^), ao = sgn(feo • e^), and 

2£. _ / j2j 



d^fci = / d^fci , (251) 

ai -/ a.ikz Gz >0 

7i"(fc)=7i'^(p), (252) 
P ^(fc|fci) = ^7' ^(p|pi). (253) 

The integration on fei is performed on the hemisphere defined by aifci • > with 
Oi = ±: 

-fci-e. >0 G [0,7r/2], 

-fci-e3>0 <^ 6li e [7r/2,7r] , ^ 
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with ki ■ = cos^i. If we sum over all the propagation directions and define the 
tensor, 

G«sv»{R,k\RQ,ko)= GTs7»iR^k\Ro,ko) (255) 

a— it ,ao— it 

which does not distinguish if the wavevectors k and fco are vertically directed upward 
or downward, the radiative transfer equation takes its classical form pi 1301 IK7j : 

k ■ \7G^sv^{R,k\Ro,ko) + KeG.^sv»{R,k\Ro,ko) 

S{k - fco) 



{Any 



■d{R-Ro)I±{k)<E>IJk) + J A^k^V{k\ki):G^sv»{RMRoM) 

(256) 

with 

7_L(fc) =1 kk, (257) 
P'\k\k,)= J2 ^'"'"'(felfci), (258) 

a— ± ,ai— ± 

^j^nw{Kk-K%)cl\K'MK'M)®cl'\K'MK'M). (259) 

Thus, the tensor G <^sv•:$>{R^ k\RQ, fep) is the Green function of the radiative transfer 
equation where the source term is discrete, localized on Rq, and emitting in the 
direction fco ^12011^' However, it is practically better to use the Green function 
^^5y-^(-R, fejilo, ^o) because we have to distinguish if the waves are upward or 
downward to write the boundary conditions. If we use the wave vectors k and /sq 
to write these boundary conditions, then 

G^^s7»{Z = 0, fc|Zo, fco; P) (260) 
= / dfci^'<5>(fc|fci) :^<iy>>(^ = 0,fci|Zo,feo;P), (261) 

gTs7»{Z = -H, fc|Zo, fco; P) (262) 
= / dkK<s>ik\ki) ■■ G\<sv°»{Z = -H,k,\Zo,ko;P) ,(263) 



with 



<<s>(fc|fei) = ^^^^<<s>(p|Pi; = 0) , (264) 



<<<,>(fe|fci) = ^^Mn'^^,^(^p\p^.p ^ 0) . (265) 

The relationship with the phenomenological radiative transfer theory is even clearer if 
we rewrite the radiative equation of section pi|l in terms of k by using the procedure 

that we have applied on G^sv-»'- 

k-yt^,{R,k)+K,t^,{R,k)^ J2 I d^fciP'"'"'(fc|fci):X^"^(i?,fci)+5^"(-R,fc),(266) 

ai = ± 

where S^(R, k) is the source term defined by: 

Sl^iR,k) = E / d'fci^'"'"(fc|fci) ■■tZARM) , (267) 



and X^^gj;(i?, fci) the reduced intensity. 
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13. Muller matrix and tensorial product 

Until now, to take into acount the polarization of the incident and scattered waves, we 
have used the tensorial product (8). However, in the optical and radar communities, 
it is customary to use Stokes vectors In introducing the product between 

two vectors and E'^ by 



E^ &E^ 



T7>1 T7'2* I 771I 7712* 
Ey Ejj + Ejj Ey 



\ 



(268) 



V {El.El*-E\,E^y*)/i J 



for E'^'^ — Elf'ev + E^f'e.n, the modified Stokes vectors of the incident and scattered 



waves are given by 



[I^\p,)]^e,E'\p,)QE'\p,)^e, 



|i?°'(Po)Hp 
2Jic{E^\p„)v E^\p^Yh) 
\ 2lui{E^\p^)v E^^p^Yh) j 

( \Efy{r)v\^ \ 

\E%{r)H? 
2Jie{E%{r)y EfyirY^) 

\ 2lui{Efy{r)y EfyivYn) j 

To transform our previous equation using the tensorial product, in equations using 
the Stokes vector, we introduce the transformation Tr*^^® by 



[I%{r)\ =e^E°^{r)QE^\r 



eo 



(269) 



(270) 



I3,I3'=H,V 



\I\ 



\ 



V i^VH - £hv)/'^ ) 

With this transformation, we easily check the following property 



Tr®^® iE" 



E^ QE^ , 



since for 



we have 



l3=H.y 



E 



2 * 



(271) 

(272) 
(273) 

(274) 



J2 ElEpep®e*p,. 

I3,I3'=H,V 

We can generalize the product and the transformation Tr®^® to dyadics operators 
/ and g 

E fl3l3'^l3e0' , g= E gpp'epep, , (275) 

13,I3'=H,V P,I3'=H,V 

acting on the vectors E^ and E'^ defined by equation H273|l . in requiring the following 
properties: 

(7 g) • {E' E^) - (7 • E') {g ■ E') , (276) 
70g = Tr®-'®(70r) , (277) 
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where the product • in the left hand side of equation H27fci|l is the usual product between 
a vector and a matrix. Then, we obtain 



fQg 



( fvvgvv 

fHVOHV 
[fvvg'HV 

+fHvgvv\ 
-'Afvvg'Hv 
\ —Invgyv] 



fvngvH 

fuHgHH 
[fvHg'HH 

+ fHHg'vH] 
-'AfvHg*HH 

— jHHgvH] 



^ Ifvv gy H~^fv H gvv] 
^ [/ff V Sff jf +/ff ff Sff v ] 

ilfvv gHH+fnH gvv 

+ fv H g'tiv + f HV g'v h\ 
iilfvvg'HH^fHHgvv 
+ fv H gHV ~ f HV g'v h] 



^[fvvgvH~fvHgvv\ ^ 

^[f HV g'u H ~ S H H g'uv^ 

i[fvvgHH~fHHgvV 

— fvHgnv+fHV gvfj] 

^IfvvgHH+fHHgvv] 

—fvHgHV~lHvgvH] / 



and 

f},P',f3o,P'a=H,V 

where 
[M] = 



(278) 
(279) 



V 



■M.HH:VV 
[J^V H;VV 
+MhV:Vv] 

-i[MvH,vv 

— J^HV:Vv] 



■Mw.HH 
■M-H H.HH 
[A4VH;HH 

-\-Mhv-hh] 

-i[MvH:HH 
-MhV;Hh] 



^[MvV-VH+MvV'Hv] 
^[MhH:VH^MhH;Hv] 
■| [J^V H:V H-\-J^HV:HV 

+ MvH;HV+MhV;Vh] 
^ [-Mv H-V H — J^HViHV 

+ -Mv H;HV — -Mjf V; V/f] 



^[Mvv-VH-Mvv-,Hv] \ 

^ H H ■,V H — J^H H:Hv] 
^ [ VH:VH — A4hV:HV 

— ■MvH;HV-\--MhV;Vh] 
■| [A4v H:V H-'rJ^HViHV 

— A4vH\HV—-M.HV;Vh] j 



(280) 

As was demonstrated in section the scattered field Ff^^ in the medium is 
expressed as a function of the incident field with the help of the operators /(p|Po): 



-fip\p,)-E^>'\p,). 



(281) 



The Stokes vector of the scattered wave can be related to the Stokes vector of the 
incident wave by the Muller matrix [M] : 

1 



[Isvir)] 



;[Mip\p,)]-[l"\p,)] 



(282) 



and from the properties in (|276|l . we show that the Muller matrix is given by the 
product © of the operators / with itself: 

[M(p|po)] - 7(p|po) 7(p|po) . (283) 

Furthermore, we have defined a tensorial cross-section which verifies 



47r — — 

'^(PoIPo) = -J /(pIpo) /(pIpo) ■ 



(284) 



With the transformation Tr®~*®, we can define a matrix cross-section [cr(p|pQ)] by 

[W{p\p,)]=Tr^--^{W{p\p,)} , (285) 
and from the properties in 1277() and equations H283I284|I we have: 



Mp\Po)] = ^ [M{p\Po) 



(286) 



Hence, the field scattered by the random medium is entirely characterized by its 
Muller matrix [iW(p|pQ)], or equivalently by [c(p|po)] which is the Muller matrix 
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representation of the tensor (t(p|pq). Moreover, we can rewrite all the equations 
previously obtained with the help of MuUer matrices in transforming the product (g) 
by the product since we have 

Tr®^® (m^ : M^^ =M^-M^ (287) 

for 

Tr®^® (ai''") = m''' . (288) 

The only non trivial transformation that we have to carry out is the right transpose 
of a tensor that we must replace by the right transpose of a MuUer matrix and which 
is defined by 

[Tr®-® {M{p\po)}f- = Tr®-® {[M{p\po)f-} . (289) 

The explicit expression of the right transpose of a MuUer matrix is given in 
( [Appendix E| |. 

14. Recapitulation 

In this section, we rewrite all the equations previously obtained in using the product 
© and the simplifications described in section H12(l . 

14-. 1. Scattering operators for the particles and the rough surfaces 

For an incident wave inside medium 1: 

E^\r) ^ E^'{ko)e''^' , (290) 

the scattered field produced by a particle of permittivity — — ei + surrounded 
by a medium of permittivity gg is described by a T matrix : 

E^%r)^ J d^rid'r2GT{r,ri) ■tl\ri,r2) ■ E\r2) , (291) 

where the Green function Gi is 

G^ir ro) = (^7+ -i^Vv) GTir - rp) , (292) 

piKe Ik-'Toll 

GTir ~ ro) = P-V.-j. ^ , (293) 

47r ||r — roll 

and K'l = K^ac- the far- field, the scattered field is a spherical wave: 

E''(.r)^^^Jik\ko), (294) 
\\r\\ 

and the operator /(A;|fco) has the following relationship with the T-matrix: 

47r7(fe|feo) = (T - kk) ■ to{K,k\Keko) ■ (7 - fcofco) • (295) 

For a spherical particle, an exact expression of /(fc|fco) is given by the Mie 
theory 
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14-2. Effective medium 

The effective permittivity Ce is solution of the foUowing system of equations: 

e,Kl^J = ei Kl^J + nCl\ka\k^) , (296) 

Cl\k\ko)=to\k\k^)+n J ^ h{k-ki)tl\k\k,)-GTiki)-Cl\k,\ko), {297) 
where 

h{k - fei) = J dV exp(-i(fe - fci) • r) [g{\\r\\) - 1] , (298) 

GTik) = J d-V exp(-ifc • r)GTir) . (299) 

c],^{k\ko) is the scattering operators of a particle located at the origin which takes into 
account the correlations with the other particles with the help of the pair distribution 
function g{r). 



14-3. Surface scattering operaiors 

The fields scattered by the upper rough surface are described by the operators 

T , R , T , and the field scattered by the lower rough surface is described by the 

operator R such that 





Incident field in the medium 


Scattered field in the medium 


r'" 













1 


r'' 


1 


1 




1 







1 


1 



For example, an upward incident plane wave in medium 1: 



E^\r) = £:i'(pQ)e'Po-="+^"-(Po)^ 

is scattered into an upward wave in medium and a downward wave in medium 1 and 
they are respectively given by 

^"'^"■^ = / (^ e'P^+'^o^P)- T°'(p|po) • E'\po) , (300) 
^''(""^ = / (^ e'^'--'«=(P)^:R°'(p|Po) • E'\po) , (301) 

where ao{p) = y/eoK^^^^ — p^, and ae(p) = \/^e^vac " P^- Since the mean plane of 
the lower rough surface is located at z = — i?, the scattering operator R^^^ describing 

21 

this surface is the scattering operator R of a surface having its mean plane located 
at 2; = multiplied by a phase factor which depends on the thickness H: 

R"^\p\Po) = e'("«(P'+"«(P''»"^:R^\p|po) • (302) 
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With these scattering operators, we can construct the scattering operators S of an 
homogeneous slab of permittivity ee having rough boundaries: 

^1-1-^—01^1+1- _i+o- ^-1+1- — 10 

s'-'^= r'' + -R'' ■r"''-R''+..., s'-'-=T'°+r-'--T'\ (303) 

-gl+l+ ^-^H21 _o+l- ^^1 -1+1- 

_0+l+^_01^^1_^l+l+_ 

For example, S^^^ described the upward wave scattered inside the slab (medium 1) 
for an incident downward wave in medium 0. 

14-4- Coherent field and coherent scattering cross-section 
For an incident field E^'^{r) in medium 

E'^\r)=E°\ko)e'^°^°-' . (304) 
The average scattered field over the random medium is specular and is given by 

« E'^{r) »sy= S'°\p^) ■ E''\k^)e^^^^-^ , (305) 

with 

Koko = - an{Po)ez , KqU = + aa{Po)ez , (306) 
< S°+°-(p|po) >5= (27r)25(p -po)5™\p) . (307) 
Then, we found for the average cross-section 

[o^'"'(fe|fco)] = 47r I cos^ol 5{k - feo)^"\po) © ^""(Po) : (308) 
where pg = -K^o feo-L and cos^q = fee • e^. 

14-. 5. Incoherent cross-section 

The incoherent Muller cross-section is decomposed in four components: 

\ = [<^L=0 J + [<^L=l \ + V^LaddeA + V^Crossed\ ■ (^09) 

14.6. Zero order scattering in volume 

The first term [ct^^q'*] described only the incoherent scattering contribution due to 
the boundaries; the random medium is replaced by an homogeneous medium with the 
permittivity Ce'- 

[^Elf{k\ko)] = 4^ co^e^l^"s> "'\k\k^) , (310) 
where €08 6* = k • Sz and S^^g^ is defined by 

{2w)H{0)Sq^s> (fe|feo) = 

x{ <S°+°-(p|po)©5"+""(p|Po) >s - <S°+°-(p|po) >s0 <S°+°"(p|po) >s }, 

(311) 

and p = Kq k±, Pq = Kq feo±. 
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14-7. First order scattering in volume 

The second component [o^"^ i'*] of the incoherent scattering cross-section contains only 
the single scattering contribution by the particles with the multiple scattering terms 
due to the boundaries: 

r / 1 ^—{a2\ scc62\—ai\ scc6i\)iie H\ 

a^^=±J'^u^2 (a2|sec6l2| - ai| sec 6*11 j^e 

47r cos -;^0+;la2 I , ;=-la2lai,,- ,f , -;^lai :0~ , , , ,„ „, 

^^'^®<^>^ '^^'^Q (fe2|fci)-50<s>(fci|fco), (312) 

where sec^? = 1/ cosf?, and = 2ImA'e. The MuUer matrix Vq^ ^ describing the 
scattering by one particle is defined by 

(4^) 



X cI\k'MK'M) cl\K'MKki) , (313) 
where Q is the Heavisidc function, and w is the structure factor of the random medium: 

wik~ko)^l + n J dV e-'('=-'=«) '- [g{r) ~ 1] . (314) 

The MuUer matrices describing the scattering at the boundaries are 

(2.)^<5(0)5°+<^-(fc|fc2) = ,/f <5°+^-(p|P2)0^""^"^(p|P2)>5, (315) 



ifg^ I cos ^1 1 —laiO-^ I laiO- 

^~(2^ 

(316) 



(27r)^^(0)5Q<s>(fci|fco) ^ a/- ' <g' (PiIPo)©'^ (PiIPo)>S, 



with 



p = Kok± , cos 6 = fc • 62 

P2=K'^k2±, cos 6*2 = A;2 ■ 62 

= K'^kij^ , cos(?i = fci • 

Po = ^ofco± , cos 00 = itQ- e.y_ 

In ( [Appendix F| |, we justify why in introducing the factors \J ^qI (-'^ a-nd ^ fe/eo in 

definitions H315|) and l|316|l . the operators Sq^^^ and <5q ^'5^ match those used in 
the phenomenological radiative transfer theory. The integral in equations IjlGSf) is 
defined by 

d2fcid2fc2= / d^fci / d^k2 (318) 

ai,a2 ^ a\h\-ez>Q J a2k2-ez>0 

where for ai — + the integration is on upper hemisphere, and for ai = — on the lower 
hemisphere: 

d'k = l /of -n^d^/g'd^ if« = +' (319) 
afc-e, >o \ /J/2 sin OdO J^"" dcf) if a = - , 
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14-8. Definition of the Green function Qq^^y^ 

To calculate the two other contributions [^Tadder] [^Crossed] 7 '^6 need the Green 
function Q<^sv^ of the radiative transfer equation. This Green function is the Wigner 
transform of the Green function G^^ describing the propagation of the electric field 



in the medium I (See paper I): 



X « gIUR + 5, iio + y) O GlyiR - y fio - y ) »5V 

- « Gly{R +^,Ro + y ) »sv « Gly{R -^,Ro- y ) »5y • (320) 
Under the quasi-uniform field approximation, we have jl3| 

K' 



Ge«sv»iR, ko) = S{\\k\\~K)S{\\ko\\~K) Gq<<sv»{R: k\Ro, feo) .(321) 



To write the boundary conditions, we have to distinguish between upward and 
downward waves, and we introduce 

^0<i;sy>>(-R>^l-Ro,fco) = 6(afc • e^)e(aofco • e^)^Q«5y»(ii, fc|fio, ^o) , (322) 

where is the Heaviside function, and a = ±, Oq — ±. Finally, as the random 
medium and the rough surfaces are statistically homogeneous, we introduce the Fourier 
transform of the Green function along the horizontal coordinates {X,Y): 

^e<<sv»iZ, k\Zo, k;P) = J d^X - X,) e-'''-(^-^°)C<<sV»(i2, k\Ro, feo) (323) 
with R = X + Zbz and Rq = Xq + ZqEz. 

14-. 9. Radiative transfer equation 

The matrices ^o<5y> are the Green functions of the radiative transfer equation: 
k ■ Vp^Q4°sy»(Z, k\Zo, fco; P) + KeGsllsv^iZ, k\Zo, ko; P) 

= ^-^J^KaJ{Z-Zo)l\\k)Ql\\k) 

+ Y. I d^kK^''\k\k,)-glZ<sv>>iz,k\Zo,ko;P), (324) 

ai=± "'"i 

where a = sgn(fc • e^), oq = sgn(feo • e^), the differential operator Vp is 

Vp-iP + e,A, (325) 
and the integral is defined by: 

d^fci = f d^fci . (326) 

'J aikz -Gz >0 
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The boundary conditions on the upper and lower surfaces are 

^0"<<sV»(^ = o,fc|Zo,fco;P) 

= / dfci^Q<5>(fc|fci).^Q<<sV>>(^ = 0,fci|2o,fco;P), (327) 

^0<sV»(-^ = -H, k\Zo, feo; P) 

dfci^Q<s>(fe|fei)-aQ<<sV>>(^ = -^^,fei|^o,feo;-P), (328) 



with 

(27r)2 5(0)<<5>(fe|fei) = <^\p\p,)Q^'\p\p,) >s, (329) 



(27r)2 5(0)<<5>(fe|fei) = ^^^^^ < r'\p\Po)QR''*{p\Po) >s , (330) 
and 

p = Klk± , cos 6 = k ■ Cz , (331) 

Po = -f^'e , COS 6*0 = feo • . (332) 

14-10. Integral equation 

The integral form of the radiative transfer equation is 

GQ'^sv>>{Z,k\Zo,ko-,P)=GT^<s>{Z,MZo,^ + ^ (4^)' / d'feid2fe2 dZsi 

xaoi"<s>(2,fc|^2i,fc2;P) •PQ''"'(fe2|fci) •e0<^^^^ (333) 
where 

gTZs>{Z, k\Zo, feo; P) = gT^\Z, k\Zo, fco; P) + Ge<s>iZ, k\Zo, fco; P) , (334) 
and 



-21 * , 



g'^t\Z,k\Zo,ko;P) = (4^)2 I eo°s^g| Q ^"."o '^a,..n(Z-Z. 



X exp (-a(Ke + ife • P){Z - Zo)/\ cos^]) , (335) 



GQ\"s>iZ,k\Zo,ko;P) = J^;^yJ^^\^T<s>ik\ko) 



X exp ^-a(K;e + ife • -P)-^/| cos6'| + ao(Ke + ifeo • -P)-^o/| cos6'o|) , (336) 

and 

(2.)^5(0)4t5>(fe|feo) = ^^^^ <S^"'"(P|PO)0^""°*(P|PO) >5, (337) 
with 

p = K'^k±, cose = k-ez, (338) 

PO = ■K'e feO-L ) COS 6o = ko- Sz. (339) 
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Ladder terms 

The third contribution to the incoherent scattering cross-section, which contains 
second order and higher scattering contributions by the particles, is 

tel.(fc|feo)]= E r dZ43dZ2i / d2fcid2fe2d2fc3d2fc4 

ai,a2,Q3,a4 = ±"'--'^ Jai, 02,03,04 

(47r)3 cos6' 7 1/11 7 1 /I n cO+'i"*/;"!/" ^ 

X— — — exp (a4 Ke^43| sec 6^4 I - Oi KeZ2l| Sect^il) 6q <g>(fe|fc4) 

I cos [74 1 

•pQ '(fc4|fc3) • C?0<<sy»(^43, fc3l^2i, fc2; -P = 0) • ' (fcalfci) • 5q <s>(fci|fco) , 

(340) 

with COS0 = k ■ e^, cos^i = fci • e^, cos04 = ^4 • and 
/ d^ki<^k2d^hd^h= I d^fei / . d^^a / . d^feg / d2fc4 . (341) 

-^01,02,03,04 J aik^-e^yO J a2kz-ez>0 Ja^k^-e^X) J a^k^-e^X) 

14-. 12. Most-crossed contributions 

The fourth contribution of the incoherent scattering cross-section, which describes the 
enhanced backscattering phenomenon, is 



[^CroLd(felfco)] = ^ / dZ43dZ21 / d^fcl d^fcj d2fc3 d2fe4 

oi,a2,a3,a4=±-^"'^ "''^i '''s ,03 ,04 

^ (4^)3 I cos go I exp(ia4Z43ife(|cos0+| - |cos04-|)-iaiZ2iif;(|cos0+| - Icos^rl)) 
I COS 6*4 I ' 

X exp ^04 ^6^43(1 sec + \ sec6^\)/2 — ai HeZ2i{\ sec0j^| + | scc9i\)/2 

S^Q<st{-^o\kA) ■ "P© ''^°'(fc4|fe3) • ^0<<SV'»(^43, fc3|2'2i, fc2; P) 



-la2lai,r- .c- , -^lai ;0- 



■Vq' ^(fc2|fci)-5Q;s>(fci|fco)J , (342) 

where the right transposed of a MuUer matrix is defined in [Appendix E| and 

P = Xo(fc± + feo±), cos0o = fco-e.. (343) 

The definition of the angles 9^ and depends on the matrix '^q^s> ^^d S^^lfgy. 
Each of these matrices has a coherent and an incoherent part: 

■^0+;loo , - I - , -^+;laoincoh f f \ , cn* t -.-c^+'Ma coh f ioaa\ 

«50<s>(fc|fco) =<50<s> (fc|fco) + '^(fc -fco)<50<s> (feo), (344) 

«50<s>(fe|fco) =<50<s> {k\ko) + S{k -ko)SQ^sy (fco) , (345) 
where A; and fco are related by the Fresnel law: 

i^ofel =i^efeoifor 5|!)t:s>(fc|fco), (346) 

i^;fcl=i^o feoffor 50<°;>(fc|feo). (347) 
For the incoherent part, the angles are defined by 



cos 0+ 1 = V 1 - (fc4± + P/2KQ^ , (348) 



cos e^l^^Jl- (fe4± - P/2KQ^ , (349) 
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cos 0+ 1 = V 1 - + P/'2K',Y , (350) 



'l-(fcii-P/2if^)2, (351) 
and for the coherent part wc have: 



cos 0+1 = VI - {Kok^/K^y , (352) 



I cos 04 I = ^1 - (Koko^/K'^)^ , (353) 

I cos 61+1 = I cos I , (354) 

Icos^jfl = Icose'+I . (355) 

We notice that in the backscattering direction, we have fe = — fcp and then P = 0, 
cos0^ = cos 6*^, cos^j*" = cos 9^. In this direction, the crossed contribution differs 
from the ladder approximation only by the presence of the right transpose around 
the MuUer matrix in the equation 1342|l . 

14-. 13. Other Formulation 

We can also calculate the incoherent cross-section in introducing the specific diffusive 

specific intensity Z'q^(Z, fc;P) which verifies the usual phenomenological radiative 
transfer equation: 

fc.VpXQ,(Z,fc;P) + ^eXQ,(Z,fc;P)= f d^kvT'^\k\k) -llZiZ, P) 

ai=± 

■•la/ 



5^'^(Z,fc;P), (356) 

""^ 



where SQ{Z,k;P) is the source term defined by 



5^«(Z,fc;P)=^ j d^k{p'^'''\k\k,)-l'^UZM-,P), (357) 



ai— ± 



and TQ^gj;(Z, fci; P) the reduced specific intensity 0: 

Jo;,rf(Zi,fci;P) =5Q<s;(fei|fco) exp(-ai(Ke+ifc-P)^i/|cos0i|) , (358) 

with cos 9i = kz -e-z- The boundary conditions necessary to solve the radiative transfer 
equation are 

t~^{Z = 0,k;P)^ [ dfei^Q<5>(fc|fei).TQ+(Z = 0,fci;P), (359) 
t^^{Z^-H,kP)^ [ dk,'R^^s>mi)-t-,{Z^-H,k:P). (360) 

Jfcie-<0 

Furthermore, we can also determine the specific intensity ll^^{Z,k; P) by using the 
integral formulation of the radiative transfer equations and its boundary conditions: 

jQ,(Z,fc;P) =^ f\z,, f A^k^d^hG'^]Zs>{ZMZ2i.kr.P)-vT''^\h\h-,P) 

AZ21 / d2fcid2fc2Ci°<5>(^,fe|^21,fe2;P) 

■ V';''"' (fcslfci; P) • (Z21, fei; P) , (361) 
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where G^q1^^s> defined by equations H334|) and H335ll33b')l . Once we have determined 
the specific intensity Iq^{Z, ki;P), the scattering cross-sections are given by 

[ol-f (fe|fco)] + rade.(fc|fco)] = 47rcos0 

X j d2fc4r|!,<5>(fc|fc4;P = o).jQ+(z = o,fc4;P = o), 

(362) 

fete<i(fe|fco)]=47r|cos0o| 

I J J cos u 

with P — k± + ko±, and {k\ko; P)] is obtained from equation H168(l in replacing 

the term 

Q—i<^2 I see 62 |— ai I see Si |)Ke -ff^ 



by 



(02! sec02| — ai| sec^iD^e 
12 (ko + 1^2 • -P)| sec(?2| — ai(Ke + ifci • P) \ sec6'i| 



(364) 



(365) 



It must be emphasized that the specific intensity Tq^{Z, k; P) is not a Stokes vector 
but a MuUer matrix. This choice was mandatory since the right transpose is defined 
only for a MuUer matrix. However, we can easily transformed the MuUcr matrix 
Xq^{Z, k; P) into a Stokes vector Xg^^^^^JyZ, k; P) with the following definition : 

^'stokes diZ,kP) -XQ,(Z,fe;P) •XQ(feo), (366) 

where XQ(fco) is the Stokes vector of the incident wave: 

X°J(fco) = E'^ko) E"^*iko) ■ (367) 

It is easy to see that the radiative transfer equation (|355|l . the boundaries condition 
H359I I36L)|) , the integral equation (|361|l and the scattering cross-section (|362|) can be 

written only in function of the Stokes vector XQ(fco)- However, this cannot be done 
for the most-crossed contribution 13631) due to the right transpose. 



15. Link with the scalar approach 



In most of the papers devoted to the enhanced backscattering phenomenon, the 
polarization and the index mismatch between the scattering medium and the 
surrounding medium is often neglected |30[ 1311 1551 f59 60, 61, 62, 69.|[70j. We can 
recover the expression obtained in this case by replacing the vectorial components 
with scalars and the product with the ordinary product between two complex 
numbers {f Q g Id*)- Without a permittiviy gradient between the scattering 
and the surrounding media (eo — — ^2) and if we suppose that the random medium 
is sparse (eg — co), we can neglect the refiexion at the boundaries, and we have 

'5©<s'>(^|fc2) = ^a2,+ '^(^ - ^2) =^ COS 6* = cos 6*2 , (368) 
^©<s>(^i|^o) = ^ai,- <5(fci - ko) =4> cos6'i = cos6'o • (369) 
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For isotropic scatterers, and in neglecting the correlations between the particles, the 
scattering operators for one particle are constant and are given, following equations 
(Ens 1213 1213) by 

t"(fc|feo) = 47r7(fe|fco) =4^/, (370) 

c"(fc|fco)~t"(fc|fco), (371) 

Vq (fealfel)- ^^-p = '^1/1 • (372) 

The effective permitivity eg and the effective constant — K^ac ^^'^ given by 
equation (|296|l : 

^ K^ + Annf , and then ^ Ki + "^fUlL . (373) 

With this effective wavevector ifg, we can obtain the extinction coefficient Kg = 
2lmKe. With these approximations and by introducing fi^ = cos 9, fii = |cos0o|, 
and the albedo of one particle a by 

a = ATTn\ff/Ke, (374) 

the single scattering contribution (|312|l by the random medium is 



incoh I 

" (375) 



Mi Ms + IM 



1 1 
\ 

Ms Mi 



where we have introduced the bistatic coefficient Yl^=i'^ which is the scattering cross- 
section cf'^£'^°^{P'si If'i) divided by the cosinus of the incident angle. Similarly, we 
derive the contribution for the ladder term from equation H34Uf) (where cos Oa = cos 9, 
cos 6*1 — cos da): 

incoh ( \ 

.incoh ( ,, \ _ " Ladder\r-s: H-i) f07a\ 
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and 





r(ri,T2,a) = (47r)a2 ^ / d% d^fca e?0<<sy>>(^43, ^31^21, ^2; P) , 



(378) 



From the integral equation H333|l on ^Q^gy~^(Z43, fc3|Z2i, fc2; -P), we obtain an 
integral equation on T(ti, T2, a): 

r(ri,T2,a) = yM^(|ri-T2|,a) + ^ j" ^ dr' VF(|ti - r'|, a) r(r', T2, a) , (379) 
where the function W is defined by 

W^(|Ti-T2|,a) = (87r) J2 I d2fc3d^fc2ao^<s>(^i,fc3|^2,fc2;P). (380) 



As we suppose that we do not have any reflections at the interfaces between the 

*0<S>1 



scattering medium and the surrounding medium, we have S Q^g-^{k\kQ) ~ 0, and we 



deduce from equation H334I336() that 



Qe:.'Ts>iZi,k\Z2,k;P)=Ge:^''\Z^,k,\Z2.k2;P). (381) 



Electromagnetic wave scattering from a random layer with rough interfaces II 44 
In neglecting the polarization in ^q^^ , the function VF(|ti — T2I, a) becomes 

^a2.sgn(Ti-r2 



W{\Tl-T2la) 

with 



E 

a2— it 



27r| COS02I 



exp 



(^-02 (^Ke 



sm O2 cos 02 
sin 02 sin 02 
cos 02 



ifc2-P (ri -T2)/Ke|cOS(?2| ,(382) 



(383) 



Since the vector P is either the null vector for ladder contribution or P = {ko±+k±)/2 
for the crossed term, then P do not have a vertical component, and P ■ k2 = P ■ k2±- 
We define angle 02 such that P ■ k2± — ||P||||fe2±|| cos02- Furthermore, we have 
sin 612 = \/l ~ cos-^ 62, and finally we obtain: 



Ifc 



2±\ 



W{\Ti-T2\,a)^ ^ exp(-Ke|ri -r2|//i2) Jo -r2|^Wil//i2^ , (384) 

where ^2 — \ cos02|, and Jq{x) is the Bessel function of order zero. 

We can also express the crossed term contribution as a function of iy(|Ti — r2|, a) 
with the help of equation (|342|l . (In this case, we have cos 9 = |cos6'^| — |cos6']~|, 



cos 6*0 — I cos 9. 
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where to obtain the last equality, we have used the fact that 7cvossed 
and r(ri, r2, a) = r*(Ti, r2, a). The factor a for the crossed contribution is given by 

llfci +fco±|| 



a 



(387) 



= [(sin^cos^ + sin0o cos6'o)^ + (sin0sin6' + sin0o sin6'o)^] /ne- (388) 

It can be easily checked that equations H375I 13771 13821 I386|l are identical to those 
obtained by Van der Mark et al. |SI1 EHl and Tsang et al. [Ml ES El- Thus, our 
approach is a generalization of the theory developed in the scalar case where the 
polarization, the boundaries, the finite size, and the correlations between the particles 
are taken into account. We could also have shown that the equations obtained by 
other groups EH |S3 EOl ED E2 EH GDI EH E] can be recovered by our approach. 
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16. Conclusion 



We have shown how to derive, starting from the Maxwell equations, the radiative 
transfer equation and its boundary conditions for a rough inhomogeneous slab. By 
using Green functions, we have obtained an integral equation, called the Bethe- 
Salpeter equation, on the intensity inside the slab. Then, by applying the Wigner 
transform to this equation and by considering the ladder approximation, we have 
obtained the integral formulation of the radiative transfer theory. The usual radiative 
transfer equation is straightforwardly obtained by differentiation. 

The main goal of this paper was to introduce the boundaries in the derivation of 
the radiative transfer equation. We have shown that this could be done by replacing 
the Green function for an infinite medium, which describes the propagation between 
two scattering events by particles, with the Green functions which take into account 
the reflexion of the waves by the boundaries. We have, in particular, demonstrated 
that the ladder contribution is identical to the phenomenological radiative transfer 
theory and satisfies boundary conditions that are identical to those derived from 
using a geometrical argument |H1 El- By starting from the Maxwell equations, 
we were also able to give an unambiguous definition of the specific intensity as a 
function of the electric field, and thus clarify the meaning of this quantity. (See 
also [HHl-) Furthermore, we were able to incorporate the correlations between the 
scatterers by multiplying the scattering operator of one particle, which appears 
in the radiative transfer equation, by a structure factor which is identical to the 
one used to describe the scattering of X-rays by crystals. Finally, we have also 
incorporated the most crossed contributions in our theory to incorporate the enhanced 
backscattering phenomenon. This contribution also satisfies a radiative transfer 
equation and boundary conditions which are slightly modified compared to the usual 
phenomenological radiative transfer theory. 



Appendix A. Dyadic and tensorial notations 

To describe the intensity of the electromagnetic field, we introduce the tensorial 
product I — E"^ ® of two vectors by 

3 

E^(g)E^ ^ El Ef, e, (g) e,, . (A.l) 

— 1 

where each vector has been decomposed on a three-dimensional orthonormal basis 
[61,62,63]: 

^1,2^^^1,2^^^ (A.2) 

4=1 

In this paper, the basis [ei, 62, 63] is either the fixed basis [e^,. By, e^], or if we know the 
wave propagation direction, the basis [e^ (p) , bh (p) , kp ] or [e\^{p), eH{p),kp ], as 
described in paper I. Then, we define the tensorial product : between two tensors 
and having the following decompositions: 

2^''' = E4''e4 0e., (A.3) 
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by 

ji:J^=^4,2S, (A.4) 

or in a equivalent way by 

{ei (g) ej/) : (e^ « e^/) = (e^ • Cfe) (e^/ • e^/) = <5i,fe(5i',fc' , (A. 5) 

where • is the usual scalar product between two vectors, and (5^ j. is the Kronecker 

I 2 

symbol. We also introduce the tensorial product between two dyads / and / having 
the following decomposition: 

by 

® f = E Mii'-jj' [eiBj (g) ei'Cj') , (A. 7) 

with 

Mu',,r ^fl,fh'- (A.8) 

If we introduce the generalized dyadic notation (e^ ® e^/ ) (e^ (g e^/ ) : 

{e-i ® e^/) (cj ® Bji) = (ciCj ® 6^/6^/) , (A. 9) 

we naturally introduce the following definition for the product between the tensor M. 
and the tensor I: 

3 

M :I= Mr,',jj'Ijf (e, ® , (A.IO) 

for 

2^ = E ^Tjy e, (g e,. , (A.ll) 

M = Y Mii'-jj' {ei ® ei>) {bj « e^/) . (A. 12) 

hi' 

This definition is equivalent to set 

[{ei (g) e^/) (Cj- (g) ey)] : (cfe (g Cfe') = (e^ (g e^/) [(e^- (g ey) : (e^ (g 6^/)] 

= (ci (g e^:)5j,k5j>^k> . (A. 13) 

Similarly, the product between two tensorial operators M and M having a 
decomposition similar to the equation ljA.12p is 

3 

AT' : Al' = ^^'^n'^3'^^o^[, (^^ ® (e,„ ® e,, ) . (A.14) 
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Appendix B. Statistical properties of the scattering operators 

It can be easily demonstrated, by using the uniqueness of the solution of the Maxwell 
equations, that the scattering operator describing an inhomogeneous slab with rough 
boundaries has the following invariance properties under horizontal translation |41| : 
^^^(P|Po)L,,,(._.),..+. = e-nP-Po)-^:R^^(p|p„)|^^^^^^^^^^ , (B.l) 

where the indexes hi^2{x — d), r-i + d mean that the rough surfaces hi{x), h2{x) and 
the particle positions ri{i — 1, . . . , N) have been translated by an horizontal vector 
d = dxe-x + dyBy. For statistical homogeneous random medium and rough boundaries, 
we have 

«^sv(p|Po)L,,,(._d),..+d »sv'=«^sv'(p|Po) ^sv, (B.2) 
p p ^ p p 

p p p p 

=« Hsy(y +Po)L, , -^-Rsvl-y +Pl - -f + Pa)\h,,,(^u, ^sv , 

(B.3) 



whatever the vector d. In using equation (|B.ip . we find 

(l_e-i(p-Po)-d) «Rsyip\p,)l^^^^^^^^^:^sv^O, (B.4) 



(^_g-i(p-p„).d) 

p p P P 

X « Rsv{^+p\^+Po)\y,,,(^)^r^®Rlv{-^+P\ - ^ +Po)L,,,,(,),,. »sy= 0. 

(B.5) 

These equalities can be satisfied, whatever the value of d, only if the average scattering 
operators are proportional to Dirac distributions: 

<^Rsv{p\PQ)h^^^(x),r, >si^= (27r)2(5(p -po)^«sv»(Po) (B-6) 
p p ^ p p 

<^Rsvi^+p\^+P0 ^ ®R*Svhl7+P\- ^+Po) , »SV 

= i2iT)HiP-Po)n^«sv»ip\Po;P)- (B.7) 
In particular, we have for P = Pa: 

— P P _* P P 

< Rsvij + Ply + Po ® Rsvi^Y + Pl - y + Po) »sy 

= (27r)2,5(0)^««sv»(p|Po;O) (B.8) 
Notice that if we suppose that the illuminated surface has a finite area A, then (5(0) 
must be defined as A/{2tt)'^. 

Appendix C. Mass operator and Dyson equation 

In paper I, we have bypassed the introduction of the mass operator and the Dyson 
equation since we have introduced from the beginning the effective permittivity in the 
definition of the Green function Gs- This mass operator is defined as 

V X Vx <G]^y{r,ro) >v -ei i^Lc <Glv{r,ro) >y= d{r - ro)! 

d3riM"(r,ri). < GsV(r-i,ro) >y . (C.l) 
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Equation (|?TT|l is referred as the Dyson equation [3 [SI |21 El ESI ■ The exact 
expression of this mass operator is the sum of all the irreducible diagrams that we 
can write and which represent the single scattering and all higher scattering process 
where recurrent scattering or correlation between scatterers are involved 56, 64 . In 
I, we have introduced the Coherent Potential Approximation (CPA) which states that 

<Gly{r,r,)>v^Gl\r,ro), (C.2) 

where G^g{r, Tq) is defined by 

V X Vx <Gl\r,ro) >v -ee i^i, < G"(r, ro) >v= S{r - ro)I . (C.3) 
In the (CPA) approach, the mass operator is, thus, given by 

M"(r,ro) = (ee-ei)'^(r-ro)7, (C.4) 
where the effective permittivity satisfied equations (177-178) in paper I. 

Appendix D. Transpose of a tensor 

In equation (|2()l|l of section Hl()|l , we carry out the right transpose of a dyadic operator 
M{p\pq) defined by 

(Ap|p4;P):P^"^^"^(p4|p3;P) 



.V \p^\p^-P):S^\s>{Px\-^V\P). (D.l) 
where P = p + Pq, and Ap = {p — Pq)/2. Since we have 

P/2 + Ap = p , -P/2 + Ap = -po , (D.2) 

P/2 - Ap = Po , -P/2 - Ap -p , (D.3) 



we deduce from the definition p69|l and H170|l of S^^l^g^ that the tensor 

M{p\pq) has the following decomposition: 

M{p\p^) = J2 Mpp.,p,p.^{el+{p) ® e°+(-po))(e;^;(po) ® e"^T{-p)) . (D.4) 

f3,P',l3o,f3l,=H,V 

From the definition of the right transpose (|187|) in the fixed basis [ex, 6^,6^] and from 
the following properties, 

4+(p) = 4-(-p), (D.5) 

enip) = -eni-p) , (D.6) 

we demonstrate that 

mp\Po)f'' = E -^J^';/3o/3i(^r(P) ® e^t(p))(e^;(po) ® e°7(po)) , (D.7) 

l3,l3',f3o,l3l, = H,V 

with 

Mfp,,p,p,=~Mp0^.^l,p, if /3V/?o, (D-9) 
for p,p\po,Po = H,V. 
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Appendix E. Transpose of Muller matrices 

To write the enhanced backscattering contribution with the help of Muller matrices, 
we need to define the right transpose of a Muller matrix. In using the transformation 
Tr®^®, we introduce the following definition of the right transpose of a Muller matrix 
[M]: 

[M]^« = Tr®^® (M^"} , (E.l) 

with 

[M] = Tr®-'® {M} , (E.2) 

and where the right transpose of the tensor Ai is defined in the previous appendix. 
Since there is a univocal relationship between a tensor M and its Muller matrix M, 
we can find the inverse Tr'^^^ of the transformation Tr®^®: 

M = Tr®^® [M] , (E.3) 

and the right transpose of a Muller matrix is 

If Miii with i,i' = 1, ... ,4 are the elements of the matrix [M], then we obtain the 
following expression for [Af]-^": 

[Mf" = 

/ Mn i [M44 - M33 i [Afi3 - M31 f [M31 + Mi3 \ 

-i M34 - i M43] -i M41 + i M14] +i M41 + i Mu] 

i[M44-M33 M22 \[M23-M32 i [M23 + M32 

+i M34 + i M43] -i M24 + i M42] -i Af24 - i Af42] 

i[Af3i-Mi3 i[Af32-M23 i[Ar33 + Ar44 5 [M34 - A/43 

+i A/42 - i A/24] -i A/41 + i A/14] -i A/12 - i A/21] +i A/12 - i A/21] 

i[M3i+A/i3 ^[A/32 + A/23 i [A/43 -A/34 5 [A/33 + A/43 

-i A/41 -i A/14] +i A/42 + i A/24] +i A/12 -i A/21] +i A/12 -i A/21] / 

Appendix F. Definition of the scattering matrix S^^^^g^ and S'Q^^g^ 

In equations (|315|l and H31()|) . of section (' (|14.7|l ). we have, respectively, introduced 

the factors (eo/fie)^^^ ^i^d (e'g/eo)^/^ in the definition of S'q^^5> and S^^gy. These 
factors insure that we recover the usual expression of the phenomenological radiative 
transfer theory. In fact, for a semi-infinite medium with plane boundaries, we have 

s"^''ip\Po) = {27:)'Sip-po)f\p,) , (F.l) 

with 

^\Po) - tviP,) 4+(Po)ey"(Po) + 4'(Po) eH{Po)eH{Po) , (F.2) 

where (pg) and (pg) ^^^^ respectively, the Fresnel coefficients in transmission for 
the polarization TM and TE. From equations (|315|l . (|248(l and H242(l . we deduce that 

5°+<\r>(fe|fco) = Sik k,)'-^T'^'-{p) (F.3) 
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where 



and 



We' cos e'o 



p^Kok±, cos9 — k-ez, (F.5) 
p = K'^ko±, cos 6*0 = fco • 62 . (F.6) 
In particular, for the polarization V V and iJ — s- iJ, we have 



7T^0+l- ^£0 cosy 01. .|2 

{p)vv-yv = ^ — \tv (p)| , (F.7) 

e' cos 9o 



T {p)hh-mh^^ -\tHiP)\ ■ (F.i 

e' cos 6^0 



The terms J^^^ {p)vv;VV and T°^^ {p)hh-hh are the usual transmission 
coefficients of intensity for the polarization (TM) and (TE) El IHI' Thus, for 

Qi 

an incident wave in medium characterized by the Stokes vector [I ], the Stokes 

Os 

vector [I ] of the transmitted wave by the plane surface is 

= 0,fc) = ^-T^^^'{p)vvyvl'^\Z = O,fco), (F.9) 

/^"(Z = 0, fc) = ^-J^^^'{p)HH,HHl^i{Z = 0, fco) . (F.IO) 

The factor eq / accounts for the spreading of the solid angle of the specific intensity 
transmitted from the medium 1 (with the permittivity e^) to medium (with 
permittivity eo) [ElEI- Hence, we see that the introduction of the term (eo/e'e)^/^ 
in definition 1)315(1 was mandatory to recover the usual expression of the transmission 
coefficients. The same reasoning can be applied to justify the introduction of the 
factor (Cg/eo)^^^ in the definition of 3^^^^. 
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